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Evaporation of Stationary Droplets 


in High Temperature Surroundings’ 


T. W. HOFFMAN? and W. H. GAUVIN® 


An experimental study of the rate of evaporation 
of stationary droplets of water, methanol, cumene, 
pentane and benzene was carried out in an electrical- 
ly-heated 9-inch stainless steel sphere. The gaseous 
environment consisted of the superheated vapor of 
the liquid under investigation and, in all cases, the 
wall and gas temperatures were the same. Drop dia- 


meters ranged from 0.4 to 1.4 mm., and temperatures 
from 100°C. to 550°C. 


The evaporation rate of stationary droplets in 
high temperature surroundings does not appear to be 
governed by the rate of heat transfer by natural 
convection, since no dependency of the Nusselt Num- 
ber on the Grashof Number was found. Attempts to 
correlate the results on the basis of Godsave’s model 
(conduction across the boundary layer controlling) 
or Spalding’s modified approach to the latter were 
equally unsuccessful. The experimental evaporation 
rates could, however, be very satisfactorily correlated 
in terms of Spalding’s Transfer Number B. 


The correlation obtained indicates that the mech- 
anism of the evaporation process is governed by the 
effect of the evolved vapors on the boundary layer 
flow. 


_ evaporation of liquid droplets in high temperature sur- 
roundings is important in an increasing number of processes 
involving the use of sprays, such as cy clone ev aporation, spray 
drying, liquid fuel combustion and in operations carried out by 
the Atomized Suspension Technique. 

The mass transfer from liquid and solid spheres has been 
extensively investigated in studies which, in the main, fall into 
two broad classes: 

(1) mass transfer from stationary liquid droplets in an 
atmosphere, the temperature of which is only moderately higher 
than that of the droplet“:?:3.4.%, 


(2) mass transfer from liquid and solid spheres in a velocity 
field, i.e. under forced convection conditions ®7:8:9:10.11,12) 

Most of these investigations were concerned with low 
transfer rates. Recently, a marked interest has developed in 
the burning of small liquid droplets .14.15.16,17,18,19,20), in which 
it has been shown that the combustion rate depends on the rate 
of evaporation of the droplets. 


\Manuscript received January 7; accepted April 27, 1960. 

2Present address: Department of Chemical Engineering, McMaster Uni- 
versity, Hamilton, Ont. 

3Department of Chemical Engineering, McGill University and Pulp and 
Paper Research Institute of Canada, Montreal, Que. 

Based on a paper presented to the C.I.C. Chemical Engineering Con- 
ference, Hamilton, Ont., November 6-11, 1959. 


The Canadian Journal of Chemical Engineering, October, 1960 


On a fait une étude expérimentale du taux d’éva- 
poration de gouttelettes stationnaires d’eau, de mé- 
thanol, de cuméne, de pentane et de benzéne dans 
une sphére d’acier inoxydable de 9” chauffée a 
Vélectricité. L’enveloppe gazeuse autour des goutte- 
lettes consistait en vapeur surchauffée du liquide 
étudié; les températures de la paroi et du gaz étaient 
les mémes dans tous les cas. Les diamétres des gouttes 
ont varié de 0.4 4 1.4 mm, et les températures de 


100°C a 550°C. 


Il ne semble pas que le taux d’évaporation de 
gouttelettes stationnaires dans un milieu a haute 
température soit régi par le taux de transport de 
chaleur par convection naturelle puisque le nombre 
de Nusselt ne dépend aucunement du nombre de 
Grashof. Il a été impossible également d’établir une 
corrélation entre les résultats, en utilisant soit le 
modéle de Godsave (ou la conduction a travers la 
couche limite est le facteur principal) ou méme 
approche modifiée proposée par Spalding. Les au- 
teurs ont obtenu cependant une bonne corrélation 
des taux d’évaporation expérimentaux en se servant 
du nombre de transport de Spalding, B. 


On obtient une corrélation indiquant que le mé- 
canisme de l’évaporation est contrélé par leffet des 
vapeurs produites sur l’écoulement 4 V’intérieur de 
la couche limite. 


Theoretical calculation of the rate of heat transfer by 
conduction to a sphere immersed in an infinite stagnant fluid 
environment predicts that the Nusselt number, (2D/k;,), should 
be equal to 2.0. When a temperature difference exists between 
the sphere and the fluid, natural convection effects arise by 
virtue of the and the 
heat transfer rate is increased accordingly. Moreover, any 
forced motion of the fluid past the sphere further increases the 


varyimg buoyancy forces around the body, 


Nusselt number above the minimum value. 


In an analysis of the convective heat transfer problem, the 
usual procedure®!.” 
the transfer of momentum, heat and mass in the boundary layer. 
These equations, known as the Prandtl boundary layer equations, 
point to the essential feature of free convection systems, in that 
the distribution of temperature and velocity are interdependent. 
In such systems, the velocity of the fluid is due entirely to the 


is to consider the equations which govern 


action of buoyancy forces arising from variations in density or 
temperature 
solution of the temperature and velocity fields impossible. It is 


The non-linearity of the equations has made direct 
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important to realize that the velocity and temperature gradients 
in the boundary layer are modified as soon as relative motion 
between the body and the surrounding fluid is established. 


Although mathematically intractable, the boundary layer 
equations have been used to derive the dynamic similarity 
relationships existing between two different systems and have 
thereby allowed a prediction of the dimensionless groups asso- 


ciated with these problems '.*, 


Ranz and Marshall'®) made an extensive study of the rate of 
heat and mass transfer from stationary liquid drops under forced 
convection, and proposed the following correlations: 


Nu = 2.0 + 0.60 ( (Pr)... (1) 


Nu’ = 2.0 + 0.60 (Re)'/? (Sc)!8.... (2) 

In his classic paper, Froessling‘* presented an approximate 

solution for the Navier-Stokes equations which gave an expres- 

sion for Nu’ of identical form, but with a coefhcient of 0.552, 

instead of 0.600. From dynamic similarity reasoning, Ranz and 

Marshall extended these correlations to include the natural con- 
vection case: 

Nau = 2.0 + 0.60 (Gr)*!/4 (Pr)'?.... 2... (3) 

= 2.0 + 0.60 (Gr)!/4(Sc)'3.... (4) 


Inherent in the derivation of Equations (3) and (4) is the assump- 
tion that material transfer did not affect the natural convection 
phenomena; these two equations have so far received little 
experimental confirmation. 

Mathers et al) assumed that the velocity normal to the 
surface was zero and neglected the inertia terms in the Navier- 
Stokes equations. They then proceeded to solve the four bound- 
ary layer equations on an analogue computer and thereby showed 
that mass transfer augmented the density gradients around a 
body. Their experimental work with spheres at very low mass 
transfer rates provided excellent experimental evidence for their 
calculations. 


Clearly, the analysis of the effect of mass transfer on the 
hydrodynamic flow pattern around an evaporating droplet is 
probably more complex than Ackermann’s simple treatment 
would suggest’). For the case of appreciable mass transfer in 
a free convection field, the velocity normal to the surface may 
be of the same magnitude as the fluid velocity in the boundary 
layer. The momentum associated with this mass will have to be 
taken into account when one attempts a solution of the momentum 


equation. This suggests that an additional boundary condition: 
dc 

Pry = D,- aty = 0..... (5) 
dy 


must be considered. 

To circumvent some of the mathematical difficulties associ- 
ated with the differential equation, velocity, concentration and 
temperature profiles i in the boundary layer have been approxi- 
mated by various polynomial e xpansions in terms of the boundary 
layer thickness'*.?*26, thus permitting integration of the equa- 
tions, ty “ of which are the integral forms presented by 
Schlichting In analyzing the simultaneous heat and mass 
transfer from a vertical wall with free convection, Spalding?” 
and Somers'® followed this procedure and considered Equation 
(5) as one of the boundary conditions. Spalding’s a priori 
solutions show that the boundary layer becomes thicker as 
greater mass transfer occurs, thereby reducing the rate of heat 
transfer from the surface. 

Spalding’? also pointed out that most analyses of convective 
heat transfer assume that: 


i-— - SS 


and are therefore applicable only to those cases where (T,, — 1) 
is smal]. This assumption may lead to an over-estimation of the 
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buoyancy forces when the temperature differences are large. In 
his own analysis, he considered the buoyancy forces to be ata 
maximum throughout the boundary layer (i.e. A7/T,, = 1), 
although this force ranges from zero at the outer edge of the 
boundary layer to its highest value at the transfer surface. 

Godsave"®, in his treatment of the rate of combustion of 
fuel drops, has shown that, at high evaporation rates, the rela- 
tively cold vapors leaving the drop can seriously reduce the heat 
transfer rate at the drop surface. This same analy sis has been 
extended or presented in a slightly different form by Goldsmith 
and Penner“®, Marshall, Ranz®, Sleicher and Churchill 
and Spalding”. Starting with the diffusion equation obtained 
from a material helenes on the diffusing component *) and assum- 
ing radial flow of vapor only (convection effects are considered 
negligible), the evaporation rate (dm/d6) can be predicted from 
the expression: 


(dm/d0) = In (1 + cpAT/X’) / (cp/4rky) (A/ra — 
where \’ = Xd — [Gr/(dm/dé)| 


bt)... ae 


and 7, is the effective radius of the so-called stagnant film 
surrounding the drop. This analysis ignores convection effects 
and attributes all transfer to conduction, much in the same way 
as Langmuir did in his pioneer work on this subject. Convee- 
tion is considered only to reduce 1. 

From this equation, the effect of evaporation rate on heat 
transfer to droplets is shown to depend on the magnitude of 
Spalding’s Transfer Number B’ = (c, 7/2’). In this analysis, 
subject to the above limitations, it is shown that the true Nusselt 
number Nuy, for pure heat transfer (that is, in the absence of 
evaporation), should be modified to yield the actual Nusselt 
number, Nu,, for the case of simultaneous heat and mass transfer, 
as follows: 


(1 /ra — 


Nu,/Nur = E(1/re — 1/7)Ke I/m) 4)... (8) 


where E = (dm/d@) c,/4k, 


By combining Equations (7) and (8): 


In(1 + B’) 


Nu, = B Nur. 


and the actual Nusselt number yields the actual rate of heat 
transfer that is available for evaporation at the droplet surface. 
Equation (9) has not been experimentally verified and the 
problem of selecting a proper temperature at which to evaluate 
the gas properties remains unsolved. 

The purpose of the present study was to investigate some of 
the problems associated with the simultaneous transfer of heat 
and mass at high temperatures. 


EXPERIMENTAL 
Equipment 

The present investigation involved measurement of the evap- 
oration rates of liquid, droplets (0.4 to 1.4-mm. in diameter) 
suspended on thin glass fibres at the centre of a heated 9-in. I.D. 
stainless steel sphere with a 3/8-in. wall. The sphere, shown in 
Figures 1 and 2 was made up of two hemispheres and was 
designed to operate at temperatures up to 1500°F. and pressures 
up to 10 atmospheres. 

The drop-injection and holder assembly, shown in Figure | 
and 2, was fitted into a 1 1/2-in. diameter port at the top of 
the sphere. The blanked-off end of a 1/2-in. tube was drilled 
to make a centered 1/16-in. diameter hole, through which the 
hypodermic needle feed system could be inserted. Guides were 
provided in the 1/2-in. tube, so that the hypodermic needle 
could be quickly inserted or removed. For atmospheric pressure 
operation, a slip seal was provided at the outermost point to 
prevent the entry of outside air. 

A fine glass filament, approximately 70 microns in diameter, 
was secured to the bottom of the 1/2-in. tube, so that its end 
was in line with the 1/16-in. diameter hole and at a distance of 
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SECTION A-A 


Figure 1—Diagrammatic sketch of the 9-in. I.D. 


\3= RADIAL TRAVERSING 
0.0005 in. D. Pt-Pt,10%Rh 


0.003 in.0 GLASS 
DROP HOLDER 








ee MICROMETER 
2 COOLING FINS 9 in. 1.0. by > in. WALL 
STAINLESS STEEL 
SPHERE 


sphere 


assembly. 


approximately 3/4-in. below it. A very small ball (diameter 
approximately 200 to 250 microns) made up of a 3/8-in. length 
of 0.0005-in. platinum, 10% rhodium wire, was fused to the 
end of the glass fibre to provide the necessary surface area to 
hold the droplet. When holding a drop, this ball was located in 
the uppermost part of the drop. No measurements of drop 
diameter were made when the drop diameter was less than 
twice the size of the ball. 

Three 0,0005-in. platinum-platinum, 10% rhodium thermo- 
couples were provided for measuring droplet and gas tempera- 
tures. The 0.0005-in. wires were fused to 0.005-in. diameter 
wire of the same material, which in turn were insulated by 
enclosing them in small glass tubes. — lead wires were 
conducted out of the sphere through a 1/8-in. diameter, thin- 
walled stainless steel tube, and this aie assembly was 
inserted into another fixed 5/16-in. tube (Figure 1). The 
thermocouples could be moved by a screwed micrometer arrange- 
ment remote from the sphere, so that the thermocouple could 
traverse the gas atmosphere surrounding the drop or be inserted 
into the drop itself. One micrometer had a large pitch thread, 
so that it could be moved quickly in and out of the drop. The 
wall temperature of the sphere was measured by three chromel- 
alumel thermocouples welded to the outside surface. 


The sphere assembly was housed in a 16-in. X 16-in. X 12-in. 
deep stainless steel box, ‘the inside surface of which was equipped 


Figure 2—Photograph showing the entire experimental 
apparatus, without the upper heater box and insulation. 
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with nichrome wire heaters. The power to each of the two 
separate electrical heater circuits in the top and bottom of the 
box was controlled by Powerstat transformers. The final 
temperature control was effected by a Minneapolis-Honeywell 
Pyrovane proportional controller. The whole assembly 
insulated with 6 inches of vermiculite. 


was 


The details of the optical bench and photography set-up are 
shown in Figure 2. It consisted of a telescope, a light-splitter 
and a variable speed Bolex motion picture camera. Shadow 
photographs of the evaporating droplets were obtained by pro- 
jecting a low power light from the side opposite to the camera. 
For the lower temperature runs, a 35-mm. camera was used, 
and in both cases the camera lens was removed. The 16-mm. 
and 35-mm. cameras gave a magnification of approximately 1.5x 
and 3x respectively. 


Procedure 

The complete apparatus was maintained at the desired temp- 
erature for about two hours prior to a run. During this time 
adjustments were made to the heater circuits, so that uniform 
temperature conditions existed in all parts of the sphere. The 
sphere was then purged of the entrapped air by injecting approxi- 
mately 75-cc. of the liquid under study. About half an hour was 


1 






Figure 3—A sequence of photographs in the evaporation of 
a water droplet showing the relative sphericity at various 
drop sizes. 
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f the drop holder. The gene of vapor in the feed tube 
caused the drop t crease slowly in size, until it separated 
Trom the necdaic ] ne aroyt sualiv ren ained on the drop holder. 
The needle was then rc ved and the camera started 1mmedi- 
ately. Some dithcult as experienced with this technique at 
the highest temperatures, in which case, after two drops had 
ailed to be 1 t er, the system was allowed to 
equilibrate before procec g aga Betore evaporation meas- 
rements were e, a p was placed on the holder and its 
temperature recorded. If this temperature was more than 0.5°C. 
below the boiling point. the system was purged again. In most 
cases, the temperature was within 0.1°C. of the boiling point, 
and at the highest temperature (~ 600°C.) any slight increase 
in temperature above the boiling point could not be detected. 
Prior to heating the apparatus, the actual magnification on 
the film was measured by rs tographing an accurately measured 


|.6500-mm. stainless steel hypodermic needle inserted through 
the in} e. The camera speed was measured by photo- 
graphing electric clock which was equipped with a sweep 





second hand 

The photographic atives were mounted on glass micro- 
scope slides and the dr opt diameter was measured on a Nikon 
Shadow graph projection microscope, with an additional magni- 
fication of 50x. The larger droplets were slightly elliptical i in 
shape (eccentricity (4b) = 1.08 max.); they were therefore 
The diameter of the sphere equal 
to the volume of this ellipsoid is given by: 


treated as prolate ellipsoids 


Wirth these small eccentricities, this diameter is essentially the 
same as that of a sphere with the same surface area as the droplet. 
4 sequence ¢ t droplet photos shown in Figure 3 indicates the 


high degree of sphericity of the droplets over the range of 
diameters measured 
Some experiments were performed with the thermocouple at 
a fixed position in the boundary layer surrounding the drop. 
was observed to increase as the droplet de- 
creased in size. Since these measurements indicated a very thin 
layer thickness (ca. three umes the drop radius), very 
large temperature gradients existed near the drop. Even with 
Yiese fine wires. therefore, considerable heat must have been 


[This temperature 


boundary 


a 
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WATER DROPLET RUNS 





conducted along the wires to the thermocouple weld. Moreover, 
since the wires extended through an unknown temperature field, 
correction of the observed temperature could not be made ; these 
observations were therefore considered as approximate oni 


Materials Investigated 


Essentially pure water was obtained by passing once-distilled 
water through an ion exchange column. Benzene, cumene, 
methyl alcohol and n- -pentane were also studied, all four liquids 
were spectroscopic grade supplied by Eastman Kodak Co.; they 
were kept under desiceant when not in use. It was found 
necessary to remove most of the dissolved gases from the liquid 
under test prior to a run, to remove the possibility of internal 
boiling during the evaporation studies. 3) has shown 
that internal boiling can result through focussing of radiation 
inside the droplet by refraction effects. 


Simpson 


The use of a one-component system alleviated some of the 
difficulties associated with choosing correct properties of the 
vapor phase surrounding the drop. The liquid density and boiling 
point data were obtained from International Critical Tables), 
The vapor properties at the higher temperatures: pre sented a 
more serious problem. Only the properties of steam have been 
measured with any precision by Keyes. The thermal con- 
ductivity data for the organic vapors were obtained by extra- 
polating the data reported by Vines“® and Vines and Bennett®? 
on a log-log plot. The specific heat data of Montgomery and 
De Vries®® and De Vries and Collins“ for benzene and meth- 
anol respectively were used in the extrapolation to the required 
higher temperatures. 


Results 

It was found that the experimental data could be represented 
by straight lines on a plot of (diameter)*-versus-time. A typical 
plot (for water at high temperatures) is shown in Fi igure 4, For 
water droplets above §00°C., a slight decrease in slope was 
observed (less than 6% from an initial drop diameter of 1.40-mm. 
to a final diameter of 0.50-mm.). At any one temperature, at 
least three measurements of evaporation rate were made at any 
given time. Some of these observations were repeated with 
different drop holders. The measurements indicated that data 
were reproducible to + 5% and in most cases complete repro- 
ducibility was obtained. 

The evaporation rate at any drop size was found from the 
slope of the D?-versus-6 plot, by means of the relationship: 


= (4/6)px [d(D)3/d0] = (4/4)pzD (d(D)?/d0].... . (11) 


(dm/d@) 






Figure 4—Measured variation of (dia- 
meter)? with time, for water droplets 
at high temperatures. 











24 26 28 W 


The Canadian Journal of Chemical Engineering, October, 1960 





The evap 
drops ars 
in Fig 
water © 
ment al 


Treatme 


The t 


can be det 


\ (dm dé 


The h 
lated to ar 
conducti\ 
conduc ted 
\n appro 
along the 
than 2% + 
conditic nn 

Since 
atmosphe 
walls, so’ 
absorbed 
of the m 
which wi 
very sma 
droplet v 
receive if 
the radia 
larger th: 
of these s 
by the dr 


Abso 
the vario 
The met! 
been disc 
drop size 


Wat 
absorpty 
ments of 
slightly | 
all temp 
clable er 
the orga 
calculate 
to show 

Fron 


where 


Combini 


Thus 
experim« 
average 

It is 
a given 


the eva 
diamete: 
number 


"Table } 
the ADI 
of Cong 

sume’ 
mm. mi 
orders p 


The Ce 





The evaporation rate data for 0.70-mm. and 1.0-mm. diameter 


drops are given in Table 1* and the results for water are plotted 
in Figure §. The data of Kobayasi"” and Nishiwaki for 
water droplets evaporating in a high temperature air environ 
ment are shown by the broken-line curves 


Treatment and Discussion of Data 


The total interfacial heat flux from the gas to the droplet 


can be determined from the following energy ‘balance 


\ (dm /d0 Ir GQ, = ky A (dT /d) = h, AAT (12) 


The heat conducted along the support-rod cannot be calcu 
lated to any degree of accuracy. Asa result of the small thermal 
conductivity and small diameter of the glass fibre, the heat 
conducted into the drop through the holder will be very small. 
An approximate calculation, similar to that used for conduction 


40 


along thermometer wells showed this heat flux to be less 
than 2°, of the total evaporative heat requirements at the worst 


condition. It was therefore ignored. 


Since in this work the droplets were suspended in a gaseous 
atmosphere which was at the same temperature as that of the 
walls, some of the thermal radiation emitted by the gas was 
absorbed in the relatively colder boundary layer. An estimate 
of the magnitude of this absorption (by means of a method 
which will be discussed in a subsequent paper) showed it to be 
very small. Therefore the thermal radiation received by the 
droplet will be essentially the same as that which it would 
receive if it were in a non-absorbing medium, 
the radiating surfaces surrounding the droplet is very much 
larger than that of the droplet itself, data on the emMiSSIVITICS 
of these surfaces are not required. Hence the radiation absorbed 
by the droplet 1s given by 


gr = 0,0A(Ty' — T,'). (13) 


Absorptivity or complete spectral transmissivity data for 
the various liquids were not available at the time of this study. 
The method for determining the absorptivity of a droplet has 
been discussed elsewhere, and can be shown to depend on 
drop size and source temperature. 


Water is decidedly opaque to thermal radiation and an 
absorptivity as high as 0.95 has been reported, The measure- 
ments of Brown) on water films 0.05-mm. thick indicated 
slightly lower values. Use of an absorptivity equal to 0.90 for 
all. temperatures and drop sizes will therefore not lead to appre- 
ciable error. A realistic absorptivity could not be assigned to 
the organic liquid droplets and therefore all these results were 
calculated twice, on the basis of a, = 0 and 1.0 respectively, 
to show the effect of this parameter 

From Equation (12) 


Nug = haD/k; = (dm/d@) (c,)/(7Dk,B’) (14 
where B’ = (c,AT)/[X — qr/(dm/d8)) . (15) 
Combining Equations (9) and (14): 

Nur = (dm/d@) (c,)/(wDk;) in (1 + B’) (16 


Thus the Nusselt numbers were calculated directly from the 
experimental data. The gas properties were evaluated at the 
average film temperature. 

It is interesting to note that, since, for a given liquid and at 
a given temperature: 


d(D)?/d@ = constant, 


the evaporative rate is directly proportional to the droplet 
diameter, as Equation (11) shows. However, the Nusselt 
number for heat transfer by natural convection is known to 


*Table 1 of this paper has been deposited as Document No. 6388 with 
the ADI Auxiliary Publications Project, Photoduplication Service, Library 
of Congress, Washington 25, D.C. A copy may be obtained by citing the 
Document No. and by remitting $1.25 for photoprints, or $1.25 for 35 
mm. microfilm. Advance payment is required. Make cheques or money 
orders payable to: Chief, Photoduplication Service, Library of Congress. 
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Figure 5—Evaporation rate data for water droplets. 


depend on (Gr)", where Ranz and Marshall suggest n = 0.25, 
as in Equation (3) while Kyte et al “ report = 0.37. One 
would therefore expect the evaporation rate under these condi- 
tions to be a function of D®” + '. However, the observations 
of this study imply that 7 = 0; this non- dependence on the 
Grashof number becomes immediately obvious if one es 
to plot the Nusselt number (true or actual) against the Grashof 
number. 


The use of [(T'x T)/T eo |* in the Grashof number 
may not satisfactorily account for the buoy ancy forces at these 
temperatures. This ratio was therefore set equal to unity 
throughout the boundary layer when the values of Gr were 
calculated, in accordance with the suggestion of Spalding”. 


\n attempt to correlate the experimental results by means of 
Spalding’s method'*! is shown in Figure 6. These plots show 
that a straight line satisfactorily represents the data for each 
liquid investigated at any chosen value of droplet absorptivity. 
Since the lines are nearly parallel, one is tempted to try to bring 
all the data together through the use of an appropriate exponent 
on the Temperature Number. However, the data at other drop 
diameters would then not correlate with these results, since the 


ordinate is a function of D-®-75, and the abscissa contains a 


radiation term which increases as /)? increases. In addition, a, 
is expected to increase with droplet diameter. The foregoing 
plots seem to indicate that the buoyancy forces are not repre- 
sented by the Grashof or Gravity number. In addition, the 
direct dependence of evaporation rate on drop diameter suggests 
that the radiant heat transfer contribution may not be separable 
from the heat transferred by convection and conduction. That 
is to say, that the overall phenomenon depends only on the 
total heat absorbed by the droplet, no matter how this heat ts 
the gravity term and radiant energy 


transferred. Accordingly, 


‘correction”’ to the latent heat of vapourization were eliminated 


*This ratio is detined as a Temperature number and was first used by Nusselt‘) 

and later by Hermann * The Grashof number was ex oressed in two parts, a 
oa ’ g : 

Gravity number G D3,p?/? and a Temperature number Te = AT/Tx 


This was suggested, so that an exponent on Je would allow for the variation 
of gas properties with temperature 
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from the Spalding correlation and (d71/d6)(c,)/(mDk;) was 
plotted against B = c, AT/X, as shown in Figure 7. All the 
data satisfactorily correlate along a single straight line of 
equation: 


[(dm/d0) (c,)/(mDky)] (Pr)-'® = 3.2 Bo-97, (17) 


The Prandtl number was included in the correlation in an 
attempt to account for the viscosity effects and hence the 
difference between the velocity and temperature profiles in the 
boundary layer of the various liquid droplets investigated. The 
exponent is consistent with that suggested by Pohlhausen’s 47) 
mathematical treatment of heat transfer phenomena. 

This equation was obtained by a least-squares analysis using 
the water droplet data only, since it was felt that the phy sical 
property data are more reliable in this case. The fact that the 
organic liquid data are slightly lower may possibly be accounted 
for as follows: 

(1) Radiation effects. At any given value of B, the radiant 
heat transfer to the water droplets 1s much higher than that for 
the organic liquid droplets. 


5.0 
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Figure 7—Correlation of all experi- 
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mental data. The experimental results Oo 
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Figure 6—Comparison of experi- 
mental results with Spalding’s correla- 
tion for free convection phenomena 
arising in combustion of liquids on a 
1%-in. spherical burner (31), 





0-6 1-0 


(2) Unreliable transport properties for the organic vapors, 
The transport properties of steam have been measured with 
high precision®®, whereas those for organic vapors had to be 
found by extrapolating low temperature data of unknown 
accuracy. 


Kobayasi’s data”) for evaporation of droplets in an air 
atmosphe re can be predicted from the correlation given in Equa- 
tion (17), which lends additional support to its validity. These 
data are shown on Figure 7 as well. Although Kobayasi’s 
experiments were conducted in an air atmosphere, the transport 
properties of the pure evolved vapors were used throughout in 
the calculations. This was suggested by the Schlieren photo- 
graphs of Nishiwaki®, which clearly show that the droplets 
are completely surrounded by a downward- flowing shell of 
evolved vapor. The one point which lies very much : above the 
correlating curve on Figure 7 refers to a benzene drop evaporat- 
ing in a 100°C. air atmosphere. At this temperature, it may be 
predicted that the concentration of air diffusing into the vapor 
boundary layer will be relatively high and the transport prop- 
erties may therefore very well be intermediate between those 
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of benzene vapor on the one hand and air on the other. 


It would be instructive to calculate the approximate magni- 
tude of the effective boundary layer thickness and how it 
changes with temperature, using a very simplified model con- 
sisting of a spherical symmetrical model. If a vapor shell of 
thickness dr and temperature T at a distance r from the centre 
of the sphere (r4<r<m) is considered, the following energy 
balance can be established: 


jar2k(dT/dr) = (dm/d0)X — qr + (dm/dO) | c,dT.. (18) 


This expression neglects radiation absorbed by the boundary 
layer since — as mentioned earlier it can be shown that 
such absorption is quite small at atmospheric pressure“. For 
water vapor the specific heat can be considered constant with 
temperature, and the thermal conductivity variation with tem- 
perature can be approximated by a straight line of the form: 


k=k, +t. (19) 


where k, = 37.5 X 1078 (cal.)/(sec.)(cm.)(°C.) and y = 
19.5 < 1078 (cal.) /(sec.) (em.) (PC.)?. 
Equation (18) can be integrated over the temperature limits 
T, to Ty to yield: 
(1/ra 1/r,) = [4rk, In (1 + ¢,AT/X’)| /c, (dm/db) + 
(4ryN’) 


- e , aaoed "7 , _ - t 7 
(cy?) (dm/d0) jin (1 + c,AT/X’) — [c,AT/(A’ + c,AT)]}... .(20) 


The values of (dm/d@) for 700-micron water drops at 
various temperatures were taken from Figure 5 and the calcu- 
lated ratios (7/ra) are shown in Table 2 for a, = 0 and a,=0.90. 

Although the values of (7%/rz) and of the effective boundary 
layer thickness (r» — 1a) calculated from Equation (20) are 
probably of the right order of magnitude, as indicated by 
Nishiwaki’s Schlieren photographs, this equation predicts 
that the boundary layer will become thinner as the temperature 
— and hence the rate of evaporation — increases, which is 
difficult to reconcile with the greater radial velocity imposed by 
the higher vapor evolution. 
classical hydrodynamics to the problem of a burning fuel droplet, 
considered as a point source of vapor in a uniform velocity 
field“, the role played by the radial vapor velocity in the 
determination of the boundary layer thickness was clearly 
implied. It must be remembered, however, that at higher 
temperatures, the addition of cold vapor and the sensible heat 
requirements of this vapor is greater, and this in turn imposes 
a relatively larger heat load at every point in the boundary 
layer. Thus the relative temperature at any point will be lower 
than that for a similar drop at lower temperatures. As a con- 
sequence of these increased downward forces, the free convection 


TABLE 2 


EFFECTIVE BOUNDARY LAYER THICKNESS 
AROUND A 700-MICRON DIAME TER Wa ATER t DROPLE:’ T 


dm/d0 X 108 a, = 0 a, = 0.90 
Tg | (Fig. 5) ae Dd aol Bis 
ee gm./sec. 








D= 700u a | ', 

| microns "/Té | microns |  "*/"4 
150 4.40 817 3:22 875 2.50 
200 9.30 806 ya) 929 | 2.65 
250 14.40 802 2.29 | 800 | 2.80 
300 20.40 | 778 y ey 960 | 2.74 
350 27.00 743 2.84 925 2.64 
400 34.40 | 712 2.03 900 ye. 
450 42.70 685 1.96 862 2.46 
500 51.20 670 1.92 833 2.30 
550 60.50 650s 1.86 813 | 2.32 
600 70.50 625 | 1.79 790 2.26 

| 





From Kumagai’s application of 


velocity will be greater and the boundary layer thickness will be 
affected accordingly. 

Since, at any given gas temperature, the rate of evaporation 
is proportional to the diameter, the vapor velocity measured at 
the drop surface is inversely prop. tional to its diameter. This 
suggests that the momentum of the vapor at the drop surface 
remains the same at all drop diameters. Equation (26) predicts 
that for any given per ie temperature, the ratio 7,/r. remains 
constant when a, = 0; this calculated ratio decreases as the 
droplet diameter ane when radiation is absorbed by the 
droplet. 

[hese results can only be suggested from the symmetrical 
model. The actual interaction of vapor momentum and cooling 
effects requires experimental observation by Schlieren or inter- 
ferometric optical techniques. 


CONCLUSIONS 


he evaporation rate of stationary droplets in high tempera- 
ture surroundings does not appear to be governed ‘by the rate 
of heat transfer by natural convection, since no depenc lency of 
the Nusselt number on the Grashof number was found. Attempts 
to correlate the results on the basis of Godsave’s model (con- 
duction across the boundary layer controlling) or Spalding’s 
approach were equally unsuccessful. 


On the other hand, the evaporation rates appear to be very 
satisfactorily predicted by the expression: 
[(dm/d8) (cy)/(aDk,)| (Pr)-!'3 = 3.2 Bo-%7,. (kay 


over the range 0.03 <B<1.0. The accuracy of this expression 
is quite encouraging, particularly in view of the uncertainties 
concerning the phy sical properties of the vapors of benzene, 
methanol, cumene and pentane. Since, in the ev aporation range 
studied, the boundary layer consists essentially of the evolved 
vapors, the values of ¢, and & for the latter at the arithmetic 
average film temperature should be used. 


It should be noted that since: 

(dm/dO)X — qr = hy AAT.............. (12 
Nug = ha D/ky = (dm/d6) (c,)/(wDR,B’)... (21) 
B’ = (c,AT)/[A — (qr)/(dm/d6)] 


and hence 
where 
an alternative expression for the proposed correlation is: 

B’ Nu, (Pr)-!? = 3.2 Bo... vdawelll 


It must be remembered that this equation can only be considered 
valid when the gas and the containing walls are at the same 
temperature. 


It is highly significant that the two forms of the correlation 
involve either the total ev aporation rate — as shown by Equation 
(17) — or the total heat transferred (convection and radiation) 

- as shown by Equation (2 rather than the individual 
transfer coefficient or Nusselt mane (Nu or Nu’) as in most 
previously published expressions. Since the mass rate of vapor 
evolution depends obviously on the evaporation rate, the infer- 
ence may be drawn that the proposed correlation relates the 

vapor flow in the boundary layer to the system conditions 
expressed by the Transfer number B. Alth« vugh the buoyancy 
effects which undoubtedly exist in the boundary layer could 
not be accounted for in the conventional way, that is, by means 
of the Grashof number, a marked downward flow of cold 
evolved vapor will be established around the droplets, which 
will in turn determine the temperature, density and velocity 
profiles in the boundary layer. ‘The magnitude of these profiles 
will be dictated by the evaporation rate, which may be thought 
of as the controlling factor, although they are obv ‘iously inter- 
dependent. 


If the proposed mechanism is correct, it might be predicted 
that the correlation should break down if high rates of evapora- 
tion are produced in cold surroundings, that is, under conditions 
where strong convective currents cannot be produced as a result 
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ot large temperature differences. Such an experimental situation 
was investigated by Diederichsen and Hall"*), who focussed the 
radiations from a carbon are onto water droplets in still air, 


through a large condenser lens. The life-time of 700-micron 
drops varied from 40 to 100 seconds, depending on the radiant 
energy level. The rate of ev aporation was found to be directly 
proportional to the square of the drop diameter, whereas a 
dependence on the first power of the diameter was found in 
the present study. 

The present experimental data might be criticized on the 
grounds that they were obtained under unsteady state conditions, 
that is with the droplet diameter continuously decreasing, and 
that the recession of the boundary might have an effect on the 
results. In his analysis of a vapor bubble liquid boundary, 
Forster“* has pointed out, however, that the relative position 
of the fluid particles is not changed by the radial motion of the 
spherical boundary, and that the heat transfer at the boundary 
should not be affected by additional convective effects. 


The experimental data and the correlating equation indicate 
that the phenomena occurring at these drop-sizes will continue 
to appear, even when the drop diameter is very small. Normally, 
natural convection effects disappear as the Grashof number or 
the size of the body becomes very smal]. Further experimental 
investigations on the ev aporation of very small droplets are 
required. It may well be that the existence of a relative velocity 
resulting from free-fall of the droplet, or turbulence effects, 
will affect the boundary layer phenomena around the droplet 
and the resulting heat transfer in a forced-flow or turbulent 
field may be less than that reported here, as other studies have 
indicated*). Moreover, such effects probably alter the ev®pora- 
tion rate of water droplets at low temperatures and very low 
Reynolds numbers, as reported by Kinzer and Gunn”. These 
aspects will also require further study. 

It is interesting to note that the group | (dm /d@) (cp) /mDk;) | 
which appears on the left-hand side of Equation (17) happens to 
be the Peclet number for the evolved vapor, being the product 
of the Reynolds number at the drop surface in a radial direction, 
and the Prandtl number. This particular Reynolds number arises 
from the momentum balance equation in the radial direction. 
This indicates that this balance cannot be neglected and empha- 
sizes once more the importance of the radial flow. 


Acknowledgements 


Assistance from the Pulp and Paper Research Institute in the form 
of an equipment grant and of financial support to T. W. Hoffman is 
gratefully acknowledged. Recognition is due also to K. Webb and A. 
Lapointe for their valuable help in the design and fabrication of the 
experimental apparatus. 


Nomenclature 
Roman Symbols 


a = major axis of an ellipsoid. 

A = surface area of droplet, cm.” 

b = minor axis of an ellipsoid. 

B = Transfer number: c,A7/A (radiation not included), 
dimensionless. 


B’ = Transfer number: c,AT [A (dr) (dm d6)\, (radia- 
tion included ), dimensionless. 

lL = heat capacity of gas, (cal.)/(gm.) (°C.). 

dc/dy = concentration gradient normal to surface, (g./cm.*)/cm. 

D = droplet diameter, cm 

D, = mass diffusion coefficient, cm.?/sec 

E = defined by equation (13) 

g = acceleration due to gravity, cm./(sec. )?. 

Gr = Grashof number: (D*gp?GAT )/(u?). 

h = heat transfer coefficient, (cal. )/(sec.) (em.*) (°C. ). 

hs = actual heat transfer coefficient at the drop surface. 

k’, = mass transfer coefficient, cm./sec. 

k, = thermal conductivity of the gas evaluated at the average 
film temperature, (cal. )/(sec.) (em.) (°C. ) 

m = mass, gm. 

m = evaporation rate, (dm/d@), gm./sec 

Nu = Nusselt number for heat transfer: hD/k., dimensionless. 

Nu, = actual Nusselt number. 

Nu; = true Nusselt number. 

Nu’ = Nusselt number for mass transfer: k,’D/D,, dimen- 
sionless. 
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Pr = Prandtl number: c,u/ky, dimensionless. 

QR = radiant heat transfer rate to drop, cal./sec. 

% = conduction heat transfer rate along drop holder, cal,/ 
sec. 

ra = radius of droplet, cm. 

rs = radius of effective spherical boundary layer, cm. 

Re = Reynolds number: DVp/y, dimensionless. 

Sc = Schmidt number: u/pD,, dimensionless. 

ie = absolute temperature, °K. 74: at droplet surface; Ty: 
at the wall; Too: temperature of the gas at infinite 
distance from the heat-receiving surface. 

Te = Temperature number: AT/To, dimensionless. 

ty = velocity of vapor normal to and at the evaporating 
surface; cm./sec 

V = stream velocity, cm./sec 


Greek Symbols 


e. = absorptivity of drop, dimensionless. 
g = volumetric coefficient of thermal expansion, °K.~'. 
Y = defined in Equation (19). 
0 = time, sec 
A = latent heat of vaporization, cal./gm. 
d’ = X — [gr/(dm/d6)], cal./gm. 
uu = absolute viscosity, gm./(cm.) (sec.). 
p = density, gm./cm.'; p,: of gas; pz: of liquid. 
o = Stefan-Boltsmann constant: 1.355 & 107-'? (cal.)/(see.) 
(cm.?) (°K. 4). 
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Pressure Drop of Single Phase Flow Through 


Raschig Ring Type Tower Packings: 
Effect of Hole Size’ 


LIANG-TSENG FAN? 


The effects of hole size of the cylindrical tower 
packing on pressure drop of single phase flow were 
investigated. It was shown that the hole size had a 
profound effect on the correlation of pressure drop 
data. This was attributed to the difference in the 
characteristics of the spaces inside and outside the 
packing units and to the interlocking of the packing 
units in the beds. 


large variety of materials has been employed for tower 
packing i in chemical and process industries. The Riaschig 
ring, which 1s a hollow cylinder whose length is usually equal 


to its diameter, is still one of the most common types of 


tower packings :??. 
Some of the major requirements of a tower packing are “ 
1. High effective contact surface area per unit volume 
of tower packing. 
High voidage and low pressure drop. 
Low weight per unit volume. 
Random packing arrangement, leading to uniform flow 
distribution. 
. Durability and low cost. 


+N 


wa 


The hole in the Raschig ring is provided to meet the first 
four requirements of tower packing mentioned above. Dura- 
bility is the main limiting factor for hole size which can be 
made in the Raschig rings. 


The pressure drops of a single phase flow through a bed 
packed with Raschig rings have been investigated frequently, 
and the results of such ‘investigations have been included in 
several of the generalized correlations of the packed bed pres- 
sure drops :*-*.%.7.8). 

In all of the previous studies, Raschig rings employed were 
mostly commercial products, and consequently, the use of var- 
ious sizes of Raschig rings caused simultaneous changes of the 
outside diameter of the packing, diameter of the hole and their 
ratios. It is, therefore, difficult to separate or distinguish the 
effects of hole size on the performance of the column packed 
with Raschig rings from the results of previous investigations. 
The purpose of this inv estigation was to study the effects of 
hole size on the pressure drop of a single phase flow. In order 
to achieve this purpose, specially manufactured clay cylinders 
with the same external dimensions but with different sizes of 
holes were employed, and all the variables of the experiment 
except the fluid flow rate were fixed. 


Packing Material and Experimental Apparatus 

The characteristics of clay packing units employed are 
tabulated in Table 1. 
seescencocecosssescsscescencsscesesccancscescescsscscencesccecescsoccocssccenes: 
iManuscript received September 24, 1959; accepted July 4, 1960. 


2Department of Chemical Engineering, Kansas State University, 
hattan, Kansas. 
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TABLE 1 


_ CHARACTE RISTICS OF Pac CKING GU: NITS 


Outside diameter, in. 


Height, in. 


Diameter of hole, 
cu. ft. 
Surface area, sq. ft. X 10? | 


Volume, 


Weight, gm. 


Packing 


Bed depth, ft. 


No. of packing per cu. ft. | 


Void fraction 


Surface area of packing, ft.? 


in, 
x< 104 





Cuar ACTERISTICS OF BE DS 





== 
Solid | Small Hole | Large Hole 
3/4 3/4 3/4 
3/4 3/4 3/4 
0 0.197 | 0.394 
1.91 1.78 1.38 
1.84 Z.ike 2.2 
10.5 | 9.6 6.8 

TABLE 2 





| Solid | Small Hole | Large Hole 

b) S9a0 2.0 2.0 
2900 3030 3430 
0.447 | 0.462 0.526 
53.3. | 64.2 | 79.3 


‘The packing units used are also illustrated in Figure 1. 


A vertical steel pipe with a 6-inch inside diameter was used 
The packed zone of the column was 2 feet long, 


as a column. 


which was held at both ends by rigid screens to the packings 
and also to prevent the expansion of the bed. 


Ring with 
small hole 
D, = 0.197 





Ring with 
large hole 
D, = 0.394 in. 


Figure 1—Packing used in experiments 
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Figure 2—Flow diagram of experiment. 
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Figure 4—Correlation of experimental data by Ergun’s 
equation. 


Precaution was taken in filling the column to assure that 
the clay packing would not be broken. This was accomplished 
by filling the column with water before putting in the packing. 
The water provided sufficient buoyancy to prevent the packing 
from being broken as well as to furnish a uniformly random- 
packed bed. The characteristics of the beds obtained are 
summarized in Table 2 

The fluid used in this experimental work was water. Its 
flow rates were measured by a flow nozzle and a mercury 
manometer. 

Another manometer was connected to the pressure taps to 
measure the pressure difference of the extreme ends of the packed 
zone. The temperature of the water at the inlet and outlet was 
constant at 18.0° C. 

Water was pumped through the flow nozzle, and the readings 
from the manometer were taken after a steady state condition 
was reached. 

Records of pressure difference were taken from the mano- 
meter for various flow rates of water. A flow diagram of the 
entire experimental set-up is shown in Figure 2. 


Results of Experiment 


Altogether 74 runs were made. The pressure drop data are 
tabulated in Table 3*. They are also plotted in Figure 3, as a 
function of flow rates. 








Oo SOLID 
4 SMALL HOLE 





O LARGE HOLE 


15 20 30 40 50 60 7080 


6, 'b/nr — ft'x10° 


Figure 3—Experimental data: pressure drop gradient vs 
water flow rate. 
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Figure 5—Interlocking of packing unit. 


Correlations of the data by Carman and Kozeny’s equation 
and Ergun’s equations are illustrated in Figure +. The results 
of these calculations are also summarized in Table 3. The 
deviations of these experimental data from the generalized 
correlations are expected’, It can definitely be observed 
from Figure + that Ergun’s correlation follows the trend of 
the data much better than Carman and Kozeny’s correlation 
does. In other words, Ergun’s correlation 1s almost parallel to 
the trend established by the experimental data. In addition, 
Figure 4 shows that the hole size has profound effect on the 
correlation of pressure drop. 


In order to simplify the analysis of the data obtained, the 
constants in Ergun’s generalized equation were slightly adjusted 
to fit the data obtained from the solid rings as follows \: 


*Table 3 of this paper has been deposited as Document No. 6389 with 
the ADI Auxiliary Publications Project, Photoduplication Service, Library 
of Congress, Washington 25, D.C. A copy may be obtained by citing the 
Document No. and by remitting $1.25 for photoprints, or $1.25 for 35 
mm. microfilm. Advance payment is required. Make cheques or money 
orders payable to: Chief, Photoduplication Service, Library of Congress. 
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- — = 1.53 + 131 — 


CL” ¢€¢ii-—e ' DG 


This equation is shown in Figure 4 by a dotted line. 


Analysis of Data, Discussion, and Conclusion 


Characteristics of bed as a function of hole size: As listed 
in Table 2, the effect of hole size to the number of packing 
of the bed was appreciable. 
The increase in the number of packing with increase of the hole 
size may be explained as a result of overlapping or interlocking 
of the packing units as illustrated in Figure 5. The consequence 
of this increase in the number of the packing units is that the 
surface area in the bed increases more than the increase of the 
area of each individual unit. The overall void space of 
the bed increases, less than the increase of the void 
due to hole space as the diameter of hole increases. The number 
of packing units per cubic foot of the bed is represented empiri- 
cally as a function of the ratio of inside diameter to outside 
diameter as 


units contained In a unit volume 


surfac: 
however, 


) tT n 


Di" . 
a - (2 
Di\’ 
2900 + 1955 
D 


Because of this increase in the number of packings, the surface 
area and porosity of the bed are expressed respectively as 
5 Ci ? Cr 
A= n,"7D2?+fr D? - D;‘ (3 
2 D ) 2D ,2 
i a th ' rie. 
€ 1 -)? ~ n C) D? D; 4) 
+ 4 ip 2 


If it were assumed thar there 1s no increase in the number of the 
packing units in the bed, total surface area in the bed and the 
of the bed would increase parabolically 


porosity as a function 


of D 


Effect of hole size on pressure drop: As already indi- 
cated in Figure 4, the effect of hole size on the friction 
factor correlation of pressure drop based on Equation (1) was 
appreciable and the experimental friction factors were always 
in the hole size increased 
In the analysis of variance performed on the 
of friction factor with respect to the hole 


less than those calculated. Increase 
this deviation 
experimental data 
size of packing ring, F-test (the variance ratio) turned out to 
be highly This indicates that the sample came from 

Therefore, it can be concluded 
ring had a definite effect on the 
Since F-test gave no clue as to how many 
differences there were, the test of Least Significant Difference 
on the data of friction factor was performed. The test showed 
of friction factors of each hole size were significantly 


. ignincant 
pe pul lation with different means 
that the hole size of packing 


friction factor 


that means 


43 fe rent 


Since increase in the hole size in the packing unit increased 
units 


this would reduce the 


the specific surface, the mean diameters of the packing 


were reduced \ccording to Kramer 


pressure drop due to the wal] effect. Calculations according to 


the Kramer's equation, however, cannot account for more than 


§&% of deviation from the krgun’s correlation 


believed thar the characteristics 
surface area with respect to the fluid flow inside 


of the Raschig ring would b« quite different from that 
In his paper 


\s already mentioned, it 1s 


of space OF 


of the outside area of the packing recent 
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indicated that a bed’s 


Greenaway “?) 
into the ring interiors, which is called the ‘ 


capacity to accept liquid 
‘inlet availability”, is 
inversely proportional to ring size. From the definition and the 
method of measurement, described in detail by Greenaway “), 
of “inlet availability”’ the interlocking of the rings would defi- 
nitely reduce the inlet av ailability of the rings over that of the 
rings with the same hole size w hen they are not interlocking, 
It may be expected that, when the hole size is small, this 
difference of characteristic is quite appreciable, but the fraction 
of the space occupied by the holes in the bed is small. The 
reverse 1s true when the hole size becomes large. Therefore, 
it is quite possible that when the hole size is larger this difference 
of characteristic appears more appreciably in an overall sense 
than when the hole size is small. 
the reasons why the deviation from friction factor correlation 
This, however, still 
does not explain the direction of deviation, which is consistently 
lower than the calculated values. 
—— 


increased as the hole size was increased. 


i 
—_— 


‘ ms | 

It can be seen that no provisions are made in Equation (1) 
to differentiate this difference of the characteristics of the spaces 
or areas inside and outside the hole in the packing units. As 
Equation (1) indicates, however, for a certain flow system at a 
fixed superficial flow rate, the porosity and the particle diameter 


are the two major factors affecting the pressure drop. Then, 


the fact that there are some differences between the effects of 


spaces or areas inside and outside the holes of the packing to 
the flow of fluid, would mean that either the effective and appar- 
ent porosities are different, or there is discrepancy between the 
effective and apparent diameters. These differences would be 
further enhanced by an inte rlocking of the particles as illustrated 
by Figure 4, to w hich the increase in packing holding c¢ capacity 
of the bed with the increase in hole size was attributed. As it 
can be observed, this interlocking of the packing units may not 
only block the passage of the fluid through the hole, but also 
hide the part of the outside surface of the packing units into 
the hole and, thus reduce the effectiveness of a part of the outside 
area. If it were considered that the interlocked particles are a 
single particle, the result would have been an increase in the 
particle volume. The effects of interlocking of the particles 
described above may never be clearly separated to allow an 
exact quantitativ ¢ analy sis because of the interaction between 
them. If it were assumed that the interlocking of the packing 
units would cause the reduction of the effective porosity, the 
result would have been a greater pressure drop than that which 
would be predicted by Equation (1). However, the results of 
the present investigation were completely reversed. The last 
two effects of the particle interlocking i.e. both the reduction 
in the effective surface area and the increase in the effective 
volume of the particle would result in the increase in the effective 
particle diameter over the apparent diameter according to the 
definitions of particle diameter, D, = 6/S, = 6/S/ y “. This in- 
crease in the effective diameter over the apparent diameter should 
result in reduction of pressure drop according to Equation (1). 


I:ffect of hole size on mass transfer: If the analogy, f © Js, 
were accepted, the results of this investigation will imply that 
the increase in mass-transfer rate due to the increases in hole 
size will not be as great as expected from the changes of apparent 
surface area, porosity and diameter. ‘This reduction of J,-factor 
may again be explained by the fact that effective surface area 
may be smaller than the apparent surface area which would be 
used in evaluating the mass-transfer coefheients. 
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Greek I 
T 


€ 
p 
Ub 


The ¢ 





Nome 


A 


( 
D, 
D, 
D 
f 


G 


iclature 


reometrical or apparent surface area of the packing units 
in one cubic ft. of the packed bed, sq. ft. 

constant defined by equation (2). 

diameter of hole, ft. 

outside diameter of cylinder, ft. 

diameter of packing unit, ft. 

friction factor as defined by Ergun, dimensionless. 

mass flow rate of fluid, Ib. /hr.-ft.? 
conversion factor, (Ib. — mass) (ft.)/(Ib. 
height of individual packing unit, ft. 

mass transfer factor, dimensionless 
depth of packed bed, ft. 

number of the solid packing in one cubic ft. of the bed. 
number of the packing unit in one cubic ft. of the bed. 
increase of the packing units over that of the solid ring, 
as defined by equation (2). 

pressure drop, lb./ft.? 

linear velocity of fluid, ft. 

geometrical or apparent surface area of a single packing 
unit, sq. ft. 

specific surface of the packing unit, ft. 
specific surface of the solid ring, ft.~" 
geometrical or apparent volume of a single packing unit, 
cu. ft. 


— force) (hr.?) 


1 


Greek Letter 


T 


rCO MM 


3.1415 

porosity of the packed beds. 
density of fluid, Ib./cu. ft. 
viscosity of fluid, Ib./(ft.) (hr.) 
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Fundamental Aspects of Solids-Gas Flow 


Part IV: The Effects of Particle Rotation, Roughness and Shape' 


L. B. TOROBIN? 


This paper examines the various complicating 
effects which rotation, surface roughness and shape 
impose on the motion and drag of single particles 
entrained in a fluid. In general, rotation appears to 
have little effect on the drag coefficient at low 
Reynolds Numbers, but influences the linearity of 
the motion and may thus play a part in the mechan- 
ism of particle entrainment in conveyed solids-gas 
systems. Negative lift has been observed in the critical 
Reynolds Number region as the result of asymmetry 
in the wake structure. Surface roughness will cause 
an early transition to a turbulent boundary layer at 
higher Re values, which will result in an appreciable 
decrease in drag coefficient, but it may increase the 
drag at low Reynolds Numbers. 


The effect of particle shape is more complex than 
has been suggested by previously published analyses. 
The available empirical correlations are discussed 
and compared. and indicate a dependency on the 
particle-to-fluid density ratio. 


i previous papers in this series '':?:*) were concerned with 
studies related to the fluid dynamics of single spheres 
moving in steady or accelerated motion through a laminar fluid. 
It is now px ssible to consider the additional complexities which 
characterize actual solids-gas flows. The present paper will be 
devoted to an analysis of the literature dealing with the effects 
of particle rotation, roughness and shape. In most cases, the 
published studies have been confined to the behavior of a single 
particle contained in a laminar fluid. The discussion will first 
deal with macroscopic observations, i.e. the measurement of 
motion and drag, and then with the microscopic empirical and 
mathematical investigation of the fluid processes. 


VI—THE EFFECT OF ROTATION 


There has been a tendency among workers in the field of 
solids-gas flow to attribute the lowered values of the drag 
coefhicient often measured in multiparticle systems to the fact 
that the particles are rotating ‘*:*’. Rotation of particles due to 
velocity gradients and collision with solid boundaries is par- 
ticularly pronounced i with non-spherical bodies. Gasterstadt “ 
reported that wheat particles would at times rotate with a periph- 
eral velocity of th e same order as that of the alr stream. How- 


ever, observations by Pasternak with freely moving spheres “” 
suggest the ratio of the peripheral velocity (hereafter referred 
to as U,) to the mean relative velocity (U.) to be of the order 


of 5 percent. Even if spherical particles have much smaller 
rotational velocities, it 1s important to know what effect they 


have on the drag coc ficient 
SESS ERR EHH EERE SEES CESSES EE SEER EERE eee 
iManuscript received April 7; accepted June 1, 1960 


2McGill University and Pulp and Paper Research Institute of Canada, 
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and W. H. GAUVIN? 


The rotation can be resolved into directions perpendicular 
and parallel to the air flow direction. The rotation with the 
axis parallel to the flow is referred to as “screw motion’’, and it 
is a far less common occurrence than the second type, which 
has been called “top spin”. The latter is brought about by 
collision with duct surfaces and particles going at different veloc- 
ities, as well as the radial velocity gradients which predominate 
in the wall region. It includes rotation about all possible axes 
perpendicular to the air flow direction. 


Top Spin Rotation 


The aerodynamics of a spinning sphere have been little inves- 
tigated. Most of the observations have been made on spheres 
mounted in air tunnels, and the Reynolds Number region in- 
volved has therefore been restricted to large values. These 
results are valuable, however, since they offer qualitative indi- 
cations of the possible effect of rotation on the drag coefficients 
at lower Reynolds Numbers. Observations at lower Reynolds 
Numbers would have to be taken with freely moving bodies 
which would be extremely difficult since relatively small angular 
velocities would cause the particle to veer considerably from 
the rectilinear path necessary for the photographic measurement 
of the rotation. 


One of the first and best instrumented measurements on 
sphere rotation was made in 1928 by Maccoll ‘*) who reported 
tests made in a 2-ft. X 2-ft. air duct on a 6-in. sphere spinning 
on an axis perpendicular to the air stream. These showed that 
at Reynolds Numbers of 6.15 X 104 to 10.70 X 104, Cp fell 
slightly from 0.52 to 0.48 as the ratio of the equatorial speed 
U, to the wind speed U, increased from 0 to 1 (see Figure 1). 
It then increased to 0.6 when the ratio reached a value of 2, and 
it remained fairly constant thereafter. 


More recently, Davies “ , Investigating the dynamics of golf 
balls, obtained results with smooth spinning spheres dropped 
into a wind-tunnel flow. Sufficient information was given to 
allow the authors of the present survey to calculate the drag 
coefficients, and these have been included in Figure 1. Although 
the drag values for the stationary case are lower than those 
given by Maccoll, they are well within the usual scatter reported 
in this Rey nolds Number region. The results, although limited, 
seem to confirm that for large Reynolds Numbers top spin will 
have little effect on the drag coefficient for the spin velocities 
expected in most solids-gas systems. 


More complete studies of the effect of top spin on drag are 
available for the simpler two-dimensional cylinder case. Experi- 
mental measurements of the effect of rotation were first obtained 
by Lafay“®"'®) and then followed by those of Reid“®, 
Betz", Prandtl», Ackeret“®, Relf and Lavender”), 
Thom “8.19, 20, 21,22,23)) and Busemann‘*4), The most recent 
study available was reported by Swanson ?®) and the drag 
coefficients he obtained are shown in Figure 2. Although the 
rotation seems to have more of an influence on the drag of the 
cylinder than it did in the sphere studies, it still required rela- 
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tively high spinning to have an appreciable effect in the sub- 
critical Reynolds Number region. In the critical region, the 
rotation caused an earlier laminar-turbulent transition mainly by 
increasing the effective relative velocity between the surface of 
the cylinder and the free stream. This could become important 
in solids-gas systems if the upper critical Reynolds Number has 
been lowered because of the relative free stream turbulence, so 
that high rotation rates would in this case cause an even further 

drop in the Reynolds Number for transition. 

The theoretical investigation of top spin rotation is again 
hampered by the presence of the non-linear terms in the Navier- 
Stokes equations which, as was previously shown“, prevent 
solutions for the non- rotating system. Lamb ‘?® and Oseen ‘?”) 
have provided solutions for stationary spinning spheres and cylin- 
ders respectively, by ignoring the inertia terms. Garstang (28) 
and later Drazin (2) using the Oseen approximation were able 
to prove that for low Rey nolds Numbers and for systems in 
which the inertia terms are not of great significance the rotation 
would have no effect on the drag of a moving sphere or cylinder. 
Analysis at higher Reynolds Numbers have been attempted by 
Bickley °, Gustafson ‘*) and Krahn ), but the complete solu- 
tion for this system cannot be expected until it has been accom- 
plished for the non- -rotating case. 


Screw Motion 


Although pure screw motion will rarely exist in solids-gas 
systems, it will probably occur as a component of the generalized 
rotation and its possible effects on the particle dynamics are 
therefore of interest. 


Figure 2—The effect of top spin rota- Cc 
tion on the drag coefficients of 
cylinders (Swanson). 
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Wieselsberger ‘*) studied the effect of rotation on fixed 
bodies spinning about an axis parallel to the air stream. These 
were performed with (a) a projectile shape 120 mm. in diameter 
and 526 mm. in length, (b) a cylinder 120 mm. in diameter and 
364 mm. in length, and (c) a disk 150 mm. in diameter ,and 
3-mm. thick. W ieselsberger was therefore able to inv estigate 
the effect of rotation on bodies in which either skin friction 
drag or form drag was predominant. He expressed his results as: 


Cp = Coo + OCp’' (U,/U,).. (1) 
Results are given for U, = 28 m./sec., and for (U,/U.) 


ratios of from 0 to 3. He found that for shape (a) there was 
little change ull (U,/U.) = 0.5, but at (U,/U.) = 2, Cp 
increased twofold. For (b) there was no increase until (U,/U,) 

1.5 and then there was a gradual increase of Cp from a value 
of Cp = 0.84 at (U,/U.) 1.5 to a value of 0.92 at (U,/U.) 
= 3. For (c) the effect of rotation was negligible. Wieselsberger 
concluded that the influence of this type of rotation was of 
importance only with long slender bodies. It diminished as the 
form drag became predominant so that for a sphere it would 
cause only a slight increase of Cp at large values of (U,/U,). 


Luthander “#) measured the drag coefficients for spheres 
under-going screw motion in the region of the critical Reynolds 
Number. Below transition, (U,/U.) ratios of up to 2 had very 
little effect on Cp, thus confirming the tendencies noted in 
Wieselsberger’s work. 
on the critical Reynolds Number only 


There was a more noteworthy effect 
at high (U,/U.) values, 





ixio® 









but it remained essentially unchanged up to values of 1.4. 
Garstang ‘**) and Drazin ‘2*) studied analytically the effects 
of screw motion for low Reynolds Number systems, and again 
found that the Stokesian drag would remain unaffected. 
Schlichting ‘**) provided an explanation of the very complex 


three-dimensional effects ir the boundary layer of rotating 
bodies of revolution in axial low for both laminar and turbulent 
boundary layers. In some cases, drag and torque were calculated, 
and these indicated increasvd drag as the rotational velocity was 
increased. The calculations could not give the entire drag 
evaluation because of the usual inability to deal with the w ake 
region. 

~ Howarth ® discussed the laminar flow in the boundary layer 
about a sphere rotatifg in a fluid at rest. He found that the fluid 
was forced from the peles towards the equator, but the mathe- 
matical solution which he presented was invalidated by the inter- 
action with the flow from the other hemisphe re which made it 
impossible to use the boundary layer equations. Howarth showed 
mathematically that for an infinite cylinder rotating in a stream, 
the drag would be unaffected by rotation. 

Burgers ‘*7) considered the effect of centrifugal force on 
boundary aie flow. The considerations were limited, however, 
to the case where the motion outside of the boundary layer is 
free from separation effects and the application of this analysis 
to the case of a sphere is difficult. 

I}lingworth ‘**) recently calculated the skin friction force for 
a body of revolution in the early stages of suddenly-initiated 
(and subsequently constant) screw motion about the axis of 
symmetry. The dependence of the separation of the boundary 
layer of a sphere on such motion was determined as a function 
of time. The equations for the longitudinal and transverse 
components of the skin friction for the sphere are given as: 


T, = 0.84625 (up,/8)!/2U, sin(x/r) [1 + 
(2.3634 + 0.484518 cos(x c/n) é&- 
} — (0.17955 + 0.007598) sin*(x/r) 


i > 
(0.48844 + 0.3 23608) cos%*x ryt €2 4 a (2) 


2wr sin(x/r) [1 + 


0.56419 (up,/6)! l 
re 32 cos(x/r)| E + 
}(0.1077 — 0.018268) sin*(x/r) + 


(8.2548 + 0.166988) cos*(x/r)} & + tic vee 


T: 


Illingworth stated that rotation will impart a centrifugal 
acceleration to the fluid elements near the surface of the body, 
which in turn produces a favorable contribution to the pressure 
gradient upstream of the maximum section of the body, and an 
adverse contribution to the gradient downstream of this section. 
The effect of increasing the spin would seem to move the posi- 
tion of separation nearer to the maximum section of the body. 
Although this would increase the total drag value, this increase 
would be insignificant for low spin rates. The effects on the 
heat and mass transfer of the rotation may be more pronounced 
since they are more directly linked to the skin friction. 


The Effect of Rotation on Linearity of Motion of a 
Particle 


Although rotation may not be an important consideration for 
the drag coefficient, the question of its being a vital part in the 
mechanism of particle entrainment into a fluid stream is still in 
doubr. It is often felt that if the velocity field about a sphere 1 1s 
not uniform, there will be a tendency for the body to migrate 
from the low velocity to the high velocity region. ‘I ‘his tendency 
has been observed by Tollert ) who found that the concentra- 
tion of spherical particles in laminar pipe flow was as much as 
25% greater at the axis as compared with the wall section. 

The difference in concentration may be due to the velocity 
gradient in the stream; it may also result from the rotation of 
the body. The friction force of the sphere rotating with top spin 
wil] cause entrainment of the fluid which adds to the flow on 
one side of the body and decreases it on the other. This phenom- 
enon is referred to as the Magnus effect, and it would tend to 
exert a lift force on the body in the direction of the increased flow. 
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This problem had been discussed as early as 1672 by 
Newton “® and later by Magnus“). Ray leigh (4) formulated 
a classical hydrody namical approach i in a paper dealing with the 
irregular flight of a tennis ball. This has been followed by the 
many studies of the rotating cylinder cited previously, but a 
comprehensive understanding of the Magnus effect is not vet at 
hand. Garstang ‘**) and Drazin °*) have also shown that the 
use of the Oseen approximation does not predict the expected lift, 


Maccoll ‘*’ pointed out that early mathematical approaches 
to this problem have visualized the spinning sphere as producing 
a swirl of air in its vicinity, modifying the stream lines so as to 
increase the pressure on one side, and decreasing i it on the other, 
By exploring the pressure variations in the vicinity of the sphere, 
Maccoll found that this is only one result of the spinning, and 
that the actual flow which he measured did not correspond to 
the mathematical treatments. He found that the velocity of the 
fluid on one side of the sphere was very low, but the velocity 
on the opposite side was only slightly above the channel speed, 
This may be interpreted as an indication of considerable energy 
dissipation i in the fluid around the sphere. The flow in the front 
of the sphere was hardly affected by the rotation. Maccoll 
established the presence of a single vortex behind the sphere 
below the axis of rotation of the body, and he found that the 
air velocity in the vortex was relatively low. The lift experi- 
enced by the sphere was expressed by ‘the following equation: 


Lift force = (1/2) Cydp, U,2.. 2... ee ee. (4) 


where C; is the coefficient of lift. A plot of C; versus the ratio 
of the equatorial speed to the wind speed shows C, to have 
initially a negative value up to a ratio of 0.5. It then rises to a 
positive value of about 0.3 for a ratio value of 2. Maccoll felt 
that there was no satisfactory explanation for this negative lift, 
but that it could be due to the flow type actually encountered. 
Similar negative lifts at low (U,/U,) ratios had been reported 
= cylinders, and this would mean that rotation would tend 

discourage rather than be the causative factor for solids 
eae Measurements taken by Fauquet (4%) on a rotating 
sphere in the Reynolds Number region of 2 X 10° to 6 X 105 
did not show a negative lift region. 


Satisfactory explanations for this phenomenon have resulted 
from the recent studies of Krahn °3), Taneda “” and Swanson (2°), 
Taneda actually photographed the trecered boundary layer flow 
of a spinning sphere and he found that in the critical Reynolds 
Number region the rotation acted to delay transition on one 
part of the “sphere and caused it to occur ‘prematurely on the 
other. This causes an asymmetry of wake structure, as observed 
by Maccoll, which in turn produces the negative lift region, i.e. 
a limited Reynolds Number range in w hich the lift is negative. 
Taneda also performed some experiments in a turbulent fluid, 
and Figure 3 shows the marked eftect of the free stream vorticity. 
If the effective turbulence is further increased by having low 
relative velocities between the particles and the stream in stream 
moving at high velocities or by the very large intensities occur- 
ring near the duct surface, it may very w ell be that the negative 
lift region would occur at even lower Rey nolds Numbers. This 
would, in this case, prevent the Magnus effect from being a 
mechanism of entrainment. 


Vil— THE EFFECT OF SURFACE ROUGHNESS ON 
THE DRAG COEFFICIENT OF A SPHERE 


Many solids-gas systems involve particles such as catalyst 
carriers whose surfaces are not smooth, and it is important to 
know to what extent this would affect their drag behavior. 

In 1916 Pannell ¢® found that the resistance to motion of a 
6-in. sphere was very much reduced by roughening its surface. 
Similar results were reported by Jacobs“ and Flachsbart “*”?. 
It was suggested that the roughness of the surface had caused 
an early transition of the boundary layer flow so that it was 
turbulent at relatively lower Reynolds Numbers. Figure 4 
shows the results of Fage and Warsap“*) on artificially ‘rough- 
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ened cylinders, and this clearly indicates a systematic decrease 
of the upper critical Reynolds Number as the size of the sand 
grains used to peers = roughening is increased. Below a 
Reynolds Number of 2 104 there does not seem to be any 
roughness effect. fade if the grain size were increased 
further, it would be more proper to consider the bluff body as 
an irregular shape. Similar experiments were performed by 
Hoerner “? with spheres, and very similar results were obtained. 


Wadell°® and Johnstone, Pigford and Chapin‘) instinct- 
ively felt that the influence of roughness on the drag coefficient 
was a minor one for subcritical Reynolds Numbers. Direct 
wind-tunnel measurement of the effects of roughness cannot be 
obtained for low particle Reynolds Numbers, but this view 
would seem to be upheld by the extrapolation of the available 
wind-tunnel data. 


Miller and McInally ‘*®), on the other hand, found that rough- 
ening of a sphere surface decreased its fall velocity at low 
Reynolds Numbers by as much as 90%, and its resistance 
behavior approached that of acube. These results were therefore 
in direct contradiction with those previously mentioned, since 
they represented an increase in resistance. 


The effect of surface irregularities may be dependent on the 
state of turbulence of the fluid. Flachsbart “7: 5) found that for 
a very low turbulence stream, he could effect a drastic reduction 
in the critical Reynolds Number by mounting a fine (0.0008-in. 
diameter) wire ring on the surface of a sphere in front of its 
equator. He found surprisingly that if the air stream was made 
very turbulent by the use of screens, the ring would no longer 
cause a reduction in drag, but, on the contrary, would increase 
the overall resistance. Very high relative intensities may very 
well cause the roughness effects to become effective at the lower 
Reynolds Numbers. There is therefore a distinct need for the 
systematic investigation of the roughness effect in combination 
with che other complicating variables. 

\ quantitative treatment of the effect of roughness is made 
difficult by the lack of an adequate system to define roughness. 
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Figure 4—The effect of surface roughness on the drag 
coefficients of cylinders (Fage and Warsap). 
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Methods similar to those employed by Nikuradse ‘* in his work 
on roughened pipe may be of value. Nemenyi‘*® suggests a 
roughness “scale” which would depend on the shape and type 
of roughness protuberances in relation to the main dimensions 
of the body. 

The pc ossibility of a minimum value of protuberance which 
would affect the drag coefficient is suggested by Michelson‘), 
He defines an ‘“‘aerodynamically smooth” surface as one whose 
irregularities are completely submerged i in the laminar boundary 
layer. He points out that since the boundary layer is very thin, 
an aerodynamically smooth surface is also smooth in the ordinary 


sense. 


VIII—THE EFFECT OF PARTICLE SHAPE 

The shape characteristics of a solid particle can exert a 
profound influence on its ability to absorb momentum from a 
moving fluid stream. It appears that the momentum transfer 
coefficients for freely moving non-spherical bodies cannot be 
obtained from fixed body wind-tunnel tests. In some cases the 
continual alterations in the body orientation can cause large 
changes in the boundary layer and wake structure. Discrepancies 
between fixed and moving drag results have also been observed 
with particles having a fixed orientation. 


The Drag Coefficient of Regular Shapes 


One of the earliest systematic studies of the fall rate of 
een? shaped bodies in a fluid was recorded by Pernolet in 
185157), This was done with accurately machined lead shapes 
of equal mass consisting of spheres, cubes, discs, plates, and 
cylinders in a vessel six meters high. The results are inconclusive, 
however, because he included the period in which the particles 
accelerated to their terminal falling velocity. The non- spherical 
particles did not fall in a straight line, and their trajectories 

varied with their shape. This means that the drag values must 

vary continuously and unless accurate measurements are taken 
at very small distance intervals, the reported coefhicients will 
represent ¢ gross averages only. 

The advent of wind-tunnel facilities at the end of the First 
World War gave rise to a great many investigations in which 
all the drag of all manner of restrained shapes were and are still 
being determined. Summaries of the results are available (8-5 
and a very thorough and intelligently _—_ survey has 
recently been published by Hoerner *°) Although these data 
cannot directly be applied to free movement, reference will be 
made to them for important pertinent information. 

Burke and Plummer ©! supported lead cubes in a pipe in an 
upward current of turbulent air by means of a thread attached to 
the pan of a balance. The thread provided only partial restraint, 
so that the results are of some interest. ¢ compared to a sphere 
of equal volume, the ratios of the resistance for cubes suspended 
from different points on their surface to that of the equivalent 
sphere were 1.7, 1.5, and 1.7 for suspension from the face, 
corner and edge, respectively. 
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Schmiedel ‘**) studied the settling velocities of thin circular 
discs moving with maximum area perpendicular to the fall. 
Subsequent workers recognized that the coefficient of resistance 
varied in a manner similar to that of a sphere, but found it 
difficult to assign a proper parameter to take the place of the 
diameter in the equations used in conjunction with the drag 

coefficients. Attempts were also made to reduce the coefficients 
to the curve obtained for spheres by the use of an appropriate 
shape factor. Towards this end Wadell*°-*®) introduced the 
concept of sphericity, which he defined as the ratio p of the 
surface area of a sphere of equal volume to that of the actual 
particle. Since the surface of an irregular body is difficult to 
estimate for irregular particles, he dev cloped another concept, 
the degree of circularity ¢ which is the ratio of the circumference 
of a circle having the same cross-sectional area as the irregular 
particle, to the actual perimeter of the latter's cross-section. 
He then proposed equations relating circularity to sphericity, 
which can be shown, however, to be only roughly approximate. 
Wadell was able to show that the data of Pernolet and of 
Schmiedel provided some evidence of the existence of a relation- 
ship between the drag coefficient and the sphericity. 

Miller and Mclnally ‘*® contributed to the accumulating drag 
data, reporting the falling velocities of different plasticene shapes 
of equal volume. The trials were performed in a water tank 
ten ‘oct j in — and six inches square. They report little or 
no tendency for well-formed spheres to spin, whereas cubes 
would side-slip and would therefore depart from a vertical path. 
Rectangular prisms and long cylinders were found to fall with 
their axes normal to their path, and they generally maintained 
a more vertical motion than did the cubes. Miller and MclInally 
obtained values for the drag coefficient, defining it by the 
equation: 

Resistance = CpL*p,U,?........ 5 ce ee 


where L = “wetted’’ diameter of the particle. The Cp values 
were consistently higher than those given for spheres and the 
data were of an exploratory nature only. 

Squires and Squires‘) studied the settling of thin discs at 
very low Reynolds Numbers, and they confirmed an earlier 
Stokesian solution given by Oberbeck‘*® that the resistance 
encountered by these shapes was given by: 


(a) 8uU.d for discs moving with the flat face perpendicular 
to the direction of the motion, and 
(b) 5.344U.d for edgewise movement in the direction of 
motion. 
‘This compares with a value of 37uU.d for a sphere. 
Squires and Squires defined the Reynolds Number for non- 
spherical bodies to be given by 


Re, = 4mp,l’,/u (6) 
where m = mean hydraulic radius of particle = (projected area 
normal to he direction of motion perimeter of the projected 


area). 

At low values of the Reynolds Number it was found that 
bodies possessing three mutually ceils planes of sym- 
metry would not change their initial orientation while falling, 
whatever it might be. This had previously been proved analyt- 
ically by Gans“ as a consequence of the linearity of the 
hydrodynamic equations of viscous flow when the inertia terms 
of the fluid are negligible. 

Heywood ‘7 -®) determined the settling velocities for four 
non-spherical shapes, for the range Re 0.01 to Re = 1000. 
He proposed the following correlating factors to account for 
the departure from spherical behaviour 

(a) the volumetric shape factor k (sometimes referred to as 

the Heywood shape factor) 


k = (a /6)d,* d (7) 

where d, = diameter of a sphere with the same volume 
as the particle, 

d. = diameter of a circle of the same area as 


the projection of the particle lying in its 
most stable position; 
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(b) the area ratio = (projected area of the particle) /(area 

of the circumscribing circle) ; 

(c) the prismoidal ratio = (mean thickness) /(maximum 

thickness). 

Heywood presented his data as a plot of Cp versus Re, 
using & as a parameter. In the Reynolds Number region below 
100, he obtained a family of parallel curves, with the highest 
value of Cp corresponding to the highest value of k. The curves 
converge and cross over at Re = 100 and then become parallel 
with the largest values of Cp corresponding to the lowest values 
of k. 

The main objection to the Heywood factor is that it is not 
specific, so that a flat and an elongated particle can have the 
same shape factor. Furthermore, k cannot be justified on hydrody- 
namical grounds because there is no reason why, for example, a 
rod and a disc having the same & should give the same danas 
istic curve. More recently ‘*) Hey wood stated that since 
minerals were mostly flat, it was the flattened value which 
should be applied. He did not feel that k could be employed 
without further modification to rod-like particles and other 
special shapes, and that further investigation of the correlating 
parameter was necessary. 

Krumbein “® investigated the fall velocity in water of seven- 
teen forms in a 6-in. diameter lucite Pipe, 6-ft. high. The 
particles were dropped in a random position and two feet were 
allowed for acceleration. The Reynolds Number region involved 
was above 500. Krumbein found that: 

(a) sphericity is not the only factor involved in correlating 
the behaviour of non- spherical shapes 

(b) it may be used as a first approximation as he was able 
to obtain a coefficient of correlation of 0.76 between 
the sphericity and the particle terminal fall velocity; 

(c) particles falling at particle Reynolds Numbers in excess 
of 500 fell with their maximum cross-section perpen- 
dicular to the direction of flow. 

Krumbein suggested an overall index of shape which he 

referred to as the “form coefficient” K, derived from the settling 
velocities of the particles concerned. It is given by: 


K = (Actual settling velocity) / 
(Settling velocity of sphere of same volume). . .(8) 


This should not be confused with the aerodynamic diameter or 
Stokes’ diameter suggested by Watson“) and Walton (72) , Fespec- 
tively. These are defined as the diameter of a sphere of the 
same density as that of the irregularly-shaped particle with a 
settling speed equal to that of the particle. 

Particular difficulties arise in the selection of the character- 
istic dimension to be used in the Reynolds Number and of the 
area term in the drag equation. Unfortunately, there is insufh- 
cient fundamental information available to make the choice 
more than intuitive. Wind-tunnel data are generally correlated 
using the downstream length of the body for the characteristic 
dimension (except for flat shapes perpendicular to the flow), 
and the projected area perpendicular to the wind direction for 
the area term. The practice favoured with freely moving bodies 
is to use dy for both the Reynolds Number and the drag 
coeflicient. 

Reynolds Number regions are conventionally classified into 
three crude groupings. A viscous or Stokesian region extends 
from zero to about 0.1, followed by an intermediate or transition 
region which extends up to about 1000, ‘The regime above 1000 
is called Newtonian or turbulent. ‘Vhe use of the latter term is 
regrettable since it often gives the impression that the particle 
finds itself in a turbulent stream, whereas the experiments are 
usually carefully confined to disturbance-free environments. 

Pettyjohn and Christiansen ‘’® studied the free fall rates of 
isometric particles in all three regions of flow. Photographic 
methods were used to measure the settling rates in liquids of 
three different viscosities. 

In the region Re < 0.5 there was no tendency for the particles 
to orient themselves in an equilibrium position. Also, the orienta- 
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tion «1 not affect the settling rates for these isometric particles. 
At Re ~ 10, tetrahedrons and cubes assumed an orientation with 
their face in a horizontal plane. At Re ~ 20, other particles 
tended to do the same, and orientation was maintained until a 
Reynolds Number region of 70 to 300 was attained, where 
particles, starting with tetrahedrons, wobbled, spun, and then 
followed a spiral vertical path. The drag coefficient increased 
abruptly with the onset of spinning, and it eventually became 
constant. It was noted that some particles started spinning at 
lower Reynolds values than did others of the same shape, and 
this produced a considerable scatter of data in the transition 
region. 

Pettyjohn and Christiansen found that for Re < 0.05, Stokes’ 
Law may be applied with an accuracy of + 2% if the correction 
factor K was employed. K was found to be given as equal to 
0.843 /n ¥/(0.065). For the region Re = 2000 to Re = 200,000 
the following equation was found to be valid within an accuracy 
of + 4% 


= V(4/3)gd, (0 — ps) / (Copy). -.--- (9) 
4.88y. 


where Cp = 5.31 Heywood’s shape factor did not 
seem to be a good criterion of the effect of shape on Cp since 
no apparent relationship resulted between k and Cp, especially 
at high Reynolds Numbers. Petty john and Christiansen con- 
cluded that ‘sphericity is a good criterion of the effect of particle 
shape on resistance. 


Very recently Chowdhury and Fritz‘) reported a study 
which seemed to repeat the work of Pettyjohn and Christiansen. 
Although they made no mention of the latter, they proposed a 
relationship for the viscous region which is in effect identical 
to that given by the earlier work. Their results for the high 
Reynolds Number region again agree closely with the Pettyjohn 
and Christiansen relationship except at the low sphericity v alues. 


Barker, a student of Christiansen, sought to extend the 
particle spectrum investigated by the latter to encompass non- 
isometric shapes. He made use of a comparatively accurate 
photoelectric technique which utilized the light reflected from a 
particle as it passed the fields of view of two photomultiplier 
cells‘, and it was found to be superior to the photographic 
methods used by previous workers, which were said to be subject 
to error when the line of fall of the particles was not vertical. 
Referring to some of Christiansen’s unpublished data, Barker 
showed that drag coefficients obtained for discs and cylinders 
having the same sphe ricities as cubes did not coincide as would 
be expected if sphericity were a good criterion. Furthermore, 
he felt that the work of Christiansen 77), Blizard ‘and 
Richards ‘"® had indicated a dependence on particle density 
other than that accounted for by the Reynolds Number. Barker 
decided, therefore, to study the effect of shape with bodies of 
varying density in liquids varying in density from 1.13 to 8.3 
g./cc. in a region which was restricted to large Reynolds 
Numbers. 


An interesting behavior of the particles in motion was noted. 
Spheres were found to wobble from side to side (1-in. maximum 
wobbling reported with a 1-in. sphere) and the deviation from 
the vertical path was random. Flat shapes fluttered but went 
straight down and only the heavier shapes spun about an axis 
in the plane of the flat surface. In one case, a particle would 
spin sometimes and flutter at others. The aver age velocity when 
fluttering occurred was found to be 5% higher than for spinning. 
Long shapes had the following ntti ristics: cylinders followed 
a vertical path except that ‘their oscillation was of a greater 
amplitude than for shorter cylinders, the square, he xagonal and 
triangular prisms went down a wide spiral trajectory. Trials 
were also attempted in a 455-ft. smoke stack, but the major 
portion of the results had to be discarded because steady state 
conditions could not be attained, 


Ihe works cited previously did not attempt to relate the 
gyrations which the particle would perform during its movement 
to the momentum transfer, although it would be expected to 








exert an important influence. Barker assumed that the tendency 
for a particle to spin was directly proportional to the applied 
force couple, and inversely proportional to the moment of inertia 
and he felt that these were bound to affect the drag. This 
reasoning was incorporated into a dimensional analysis which 
produced the relationship: 


Cos = (p/p.)" (dmin/dmax)™ ® Re ere 


where 7, m and ® are experimentally determined. 


Barker initially employed the actual projected area of the 
solid as it fell, and the length of the body in the direction of 
motion, as his values for the area and diameter in the calculation 
of the drag coefficient. He found that the curve for y = 0.67 
crossed over the curve for W 0.806 at Re = 2,700 in the 
case of lucite, and at Re = 6,600 for the case of steel, although 
the intersection could be expected to occur at the same Re value 
if the Reynolds Number was properly defined. He then employed 
the length of the cylinder instead of the longest projected length 
in the direction of fall, and then the intersections did occur at 
the same Reynolds Number and this method was considered 
best for the purpose of correlation. The area term A, used for 
the calculation of Cp was the projected area of an equivalent 
sphere for the isometric particles, and the projected area of the 
particle perpendicular to flow for the non-isometric particles. 

The results obtained for a wide variety of shapes are plotted 
in Figures 5, 6 and 7. The low scatter seems to confirm Barker’s 
method of correlation for this Reynolds Number region. Figure 
8 clearly shows the previously -overlooked influence of the dens- 
ity ratio on the drag coeflicient. This complication is probably 
of little i importance in the intermediate and low Reynolds Number 
regions where settling has been observed to be steady. 


The drag results for non-isometric particles at low and 
intermediate Reynolds Numbers are of special interest since 
they are often of direct use in practical systems in which the 
fluid is laminar. MeNown and Malaika ‘8-8? systematically 
studied the settling of non-isometric particles in these regions 
with axisymmetrical particles which were symmetrical about 
three mutually perpendicular planes. These included circular 
cylinders, double cones, square prisms and spheroids. Here 
on the particles were stable in any settling position below 
Re = 0.1 and in the region 0.1 < Re < 500 they settled with 
the largest cross-section in the three mutually perpendicular 
planes oriented perpendicular to the motion. The runs in the 
low Reynolds Number region were initiated with the particles 
oriented in their manner of settling for the intermediate runs. 
McNown and Malaika also evaluated the Oberbeck integrals 
for spheroids settling along a principal axis, and found close 
agreement between their experimental results at low Reynolds 
Numbers and the theoretically-derived ones by approximating 
the actual shapes with their spheroid counterpart. For the 
Stokesian regime, K (the ratio of the actual settling velocity to 
the velocity of a sphere of equal volume) was closely approx- 
imated by the relationship: 

a”B” 


11) 
(0.811 /¢!5%) + 0.1898” . 


For the entire range investigated, it was found that particles 
with the same a/s ratios had the same drag coefhcient. (a is 
the particle length perpendicular to the largest circle or square 
cross-section and s is the diameter ot the maximum cross-section 
or side of square base of a square prism). This indicates that 
sphericity cannot be used as an aerodynamic criterion with non- 
isometric particles for all regimes of flow. 


Pramanik ‘* extended the investigation to deal with 
triaxial particles, i.e. a, b, and ¢, the smallest, termediate and 
largest semiaxes of the ellipsoids approximating the particles 
respectively are all different (in the previous study 6 = ¢ = 5). 
Ihe results in the Stokesian range showed very good agreement 
with Equation 11. Ifa is the semiaxis in the direction of motion, 
it was found that the variation of K with the Reynolds Number 
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Figure 6—Drag coefficients of non-spherical particles (Barker). 
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Figure 11—K,/K, as a function of Re (Presler). 


was almost a unique function of a Vbc, with the b/c ratio being 
a noticeable but secondary effect. 

Presler “using the experimental data of Malaika and 
Pramanik has dev eloped a convenient generalized method for 
calculating particle drag coefficients for the Stokesian and inter- 
mediate range. Assuming that K, can be determined by assuming 
the particle resistance is equal to the resistance of an ellipsoid 
of equal dimensions, he derived it theoretically for ellipses of 
full range of semiaxis values by using the Oberbeck integrals. 
K, is obtained from the relation: 


K, = M; + Ma (b/e — 1)......... 82) 


where M, and My are analytical coefficients. 
Vf, and Mz are given graphically from the solutions as 


functions of the ratio a/Wbc and are shown in Fi igures 9 and 10. 
For the intermediate Reynolds Number region, Presler found 
it convenient to use an overall Stokes Law ‘velocity correction 
factor K,. He determined the ratios of K, to K, from the experi- 
mental data and these are given in Figure 11 as a function of 


dpU.po 


UL 
The fact that the K, values were calculated from the 
Oberbeck integrals by assuming the maximum particle area 
being perpendicular to the flow does not affect K,, since they 
will settle with that orientation in the intermediate region. This 
is not necessarily the case in the Stokesian region, and the effect 
of orientation at very low Reynolds Numbers was investigated 
by Heiss and Coull °° with particles which had one of the three 
mutually perpendicular axes equal in length. Their results com- 
pared favorably with solutions to the Oberbeck integrals given by 
Gans ‘®) for ellipsoids of revolution (i.e. spheroids since b = ¢). 
For orientations where the circularity of the projected area 
varies with different shape, their results were approximated by: 


d,\'"? (d, 


~().25 
vi d,, 


1 
logio( K,) 1 }+logio ) we as Cee) 
dy 


and in orientations where the circularity is constant and remains 
unaltered by changing shape: 


' 0.27 d, d,\¥'? 
login K,) = - 1 + logio { - ~e G F4) 
d, 0.345 d,, 1, 
( y )i/2 
d, 
where dn = diameter of circle of area equal to projected area, 


Ludwig ‘6) investigated the fall of discs and cylinders and also 


found a marked dependency on a Vbc, but the simple relation- 
ships which he obtained are specific to his limited data. 
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Figure 12—The effect of edge round- 
10> ness and orientation on the drag co- 
efficients of square cylinders 


(Polhamus). 


Becker ‘*”) has recently proposed an overall solution to the 
shape problem by developing what he terms a quadric resistance 
“law”. In an analogue to Budryk’s unsuccessful pseudo- 
theoretical treatment for spheres (as discussed in Reference 1), 
he feels that the total drag force on a particle, R, is made up 
of a viscous term Ry and an inertia term R; (these are not to 
be confused with the skin friction drag and the form drag 

respectively). 1.e. 
ae ee (15) 


Becker feels that this is a consequence of the Navier-Stokes 
equations containing both viscous and inertia terms, so that the 
solutions should contain additive viscous and inertia terms. At 
low Reynolds Numbers the inertia terms (which would be 
represented by R,) are to be neglected and at high values the 
viscous terms, represented by Ry, are held to be ineffectual. 
Presumably, the intermediate region is found by adding the two 
terms. Although Equation (15) may have some use as a corre- 
lating device, the implied connotation is conceptually incorrect. 
Although it is true that the overall drag at a high Reynolds 
Number is proportional to the square of the velocity, it is 
incorrect to say that viscous forces play no part in determining 
the drag. Were it not for viscosity, then the inviscid flow 
equations would be valid, and these would of course predict 
zero drag. In the intermediate zone, the equation governing the 
drag at very low Re values would be altered by the inertial 
effects. For the particle “viscous drag force”, Becker uses 
the equation: 


ie : 
Ry = g CvhUd, (16) 


where the coefficient Cy) 1s a function of the form sphericity 
@ = md,*/4+A,, A, being the area of the particle perpendicular 
to the motion. Cy is evaluated from the results at very low Re 
but is supposedly independent of the Revnolds Number. In the 
intermediate region, R, is obtained by subtracting the value of 
Ry calculated from Equation (16) from the overall experi- 
mentally determined one. An inertial drag coefficient, C;, is 
defined by the equation: 


Ry = (1/2) CyA,p.U? (17) 


where A, is the area projected by a sphere with the same surface 
area as the shape. Plots of C; versus Re for the region 107 < 
Re < 10 for isometric bodies show a large scatter (as compared 
to the good correlation for all shapes cited previously). Even 
more important, they point out the theoretical weakness of 
Becker's argument, since the inertia coefficients lie below the 
curve for spheres, whereas the effect of shape — particularly 
when the particles fall with their maximum cross-section per- 
pendicular to the flow — must be to increase the inertial effects 
in the flow about the body. 

At Reynolds Numbers greater than 1,000, C; becomes 
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equivalent to Cpy, if the area term is the projected area of a 
sphere of equal volume. He suggests that C; can be calculated 
from the sphericity Y by the equation 


Cy = 2,25 + 0.283 (1 ~ OF)..650..5-4505- (18) 


Becker seems to have overlooked the works of Barker, Presler 
and Pramanik. The authors of the present survey have compared 
their experimental data with those predicted by Becker’s corre- 
lations and have found the latter to be inappropriate for non- 
isometric bodies, particularly if they are prolate or oblate. 


Particle Edge Effects 


The edges of the particles used in most of the free-fall 
experiments were usually sharp, whereas the solids in multi- 
particle systems would often have their edges smoothed by 
attrition. There have been wind-tunnel studies made of the 
effect of edge roundness for large, but subcritical Reynolds 
Numbers and these may serve as a useful indication of what 
occurs in the free particle case. Polhamus“* has recently 
studied the combined effect of body orientation and edge round- 
ness with a two-dimensional square cylinder, and his results, 
which were presented together with pertinent data given by 
Wieselsberger **’ and Delany and Sorensen‘*®) are shown in 
Figure 12. It would appear that the edge roundness effect is 
most pronounced when the body is moving broadside-on, and 
in this case the drag is gradually reduced to that of a sphere, as 
the radius of edge roundness r’ is increased. When the orienta- 
tion angle a’ approaches 45°, the effect of edge roundness 
becomes negligible. These results would have to be extended 
to three-dimensional shapes of different cross-section before any 
generalized conclusions can be drawn. Since the drag calculated 
for ellipsoids closely approximates the results obtained with 
their sharp-edged counterparts, it would seem that edge condi- 
tions would not affect the drag in the Stokesian region. This 
view is substantiated by Barth ’® who noted that the resistance 
of sharp-edged mate rials increased faster with i increasing velocity 
than did that of rounded or spherical particles. 


The Coefficient of Drag of Irregular Shapes 

The correlation of the drag coefficients for irregular particles 
is difficult because of the variety of random geometric structures 
and surface characteristics involved. The work reported is 
therefore specific to the conditions and types of material studied, 
and only the most general qualitative behavior can be described. 

Schulz “) determined the constant settling velocity of marble, 
pyrite, sphalerite, quartz and coal particles, defining their shapes 
as roundish, longish and flat. The flat bodies have slightly 
higher Cp values than that given by the standard curve in the 
Reynolds Number region of 10 to 100. The data for cubical 
materials are below the standard curve at low values of the 
Reynolds Number, and then cross over at a value of about 60. 

“The resistance of coal particles was studied by Feret ??, 
Decaux ‘“®), Audibert *, Hirst ©, and Needham “8, These 
values vary considerably, and the results seem to depend 0 on 
whether the particles were balanced in an ascending fluid, « 
were descending through a fluid at rest. 

Blizard’*) was able to correlate data for the fall rates of 
powdered coal and coke particles in hydrogen and air at Reynolds 
Numbers in excess of 1,000 by the following expression: 


U, = (d,°-*)/(3.44/p,°*") (p/p,)* (19) 


This expression indicates 
work of Barker 

Burke and Plummer “*’ suspended coke particles in a turbu- 
lent air stream and obtained the relationship 


a dependence on density as did the 


U, = 47.8 [d,(p — p.)/p.)°* (20) 


where da particle mean diameter. ‘The relationship obtained 
from the data of Martin for air velocities required to support 
of crushed quartz grains is identical in form, but 
Both expressions are for Reynok Is 


Various sizes 
with a coefhaent of 50.9 
Numbers greater than 500 
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TABLE 1 
DRAG COEFFICIENTS FOR V ‘ARIOU Ss SHAPE s 


| 
ss Coal, Anthracite, 


Re Spheres Sandstone and Py rites Shale 
10 is | 3.6 3.6 
100 | 0.4 1.2 1.6 
1,000 [ 0.35 hin | 
10,000 0.35 | 1:3 | 1.8 


Miller and MclInally ‘**) found that the data for coal, anthra- 
cite, sandstone and pyrite particles all fell (within the experi- 
mental error) on the same curve. The values for shale followed 
the same curve up to a Reynolds Number of about 100 and 
then Cp became larger for the flat shale particles than it did 
for the rounded-off or cubical particles, and it seemed to increase 
slightly with increased Reynolds Numbers, as shown in Table 1. 

The data of Prockat and Linsel® and also Rosin and 
Kaiser for various types of coal, show the Cp values for 
these irregular particles to be lower than that given for spheres, 
for Reynolds Numbers less than 100, in contradiction to the 
above results. 


The proper definition of the diameter of an irregular particle 
which should be used in the definition of the Reynolds group in 
the drag coefficient relationship is difficult to establish. Eirich (%) 
has shown that porous or indented spheres carry along some of 
the fluid through which they pass. This led him to believe that 
the outer outline of an irregular particle is the one which has 
the most hydrodynamic significance, because the retained fluid 
acts as part of the particle. Haw ksley °? feels that the drag 
of non-spherical and rough particles is governed by the size of 
the enveloping surface and he suggests as a convenient approxi- 
mation the use of a drag diameter which would be the diameter 
of a sphere having the same external surface area as the particle. 
He feels that the drag diameter would be the more fundamental 
quantity when compared with Stokes’ diameter discussed |earlier. 

Guruswamy and Roy “°*) have defined the Reynolds Number 
of irregularly shaped particles as: 


Re = (V"8U,p,)/(m)......... 


where V is the volume of the irregular particle. They were 
able to correlate the settling rates of irregular coal particles 
with a plot of K versus y. 

Brownell, Dombrowski and Dickey ') have prepared a plot 
of sphericity versus bulk porosity to facilitate the measurement 
of the effective particle area and volume for complex particle 
shapes. Nichols and Gauvin“®®, however, have shown this 
plot to be very inaccurate. 

The problem of classifying irregular particles with an 
aerodynamically meaningful shape parameter must be overcome 
before any order can be brought into the determination of the 
drag coefficients of naturally-occurring particles. Sedimentary 
petrologists are immediately concerned with shape characteriza- 
tion since it is closely related to many geological phenomena, 
and pertinent information from this field of specialization has 
been discussed by Marchildon°®, Wadell°” had, in addition 
to his sphericity and circularity concepts, defined several others 
which may be of some use. He described a laborious but 
significant method of classifying the “roundness” of a particle 
by reproducing it as a camera lucida drawing and measuring 
the various radii of curvature of the plane cross sections. The 
roundness would be quantitatively obtained from maximum 
plane cross-section by measuring the average corner radius and 
dividing it by the radius of the largest inscribable circle. Wadell 
also defined an operational sphericity ¥.» to be used in approxi- 
mating the difficulty measured y, and it is given by: 


Particle Volume yo 


Yor = - (22) 
“4 Ce of smallest circ umscribing sphere 
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TABLE 2 
3 - ~ ——— 
b/c | a/b Shape Classification 
_| ‘a Seana 
» Zs > 2/3 Sphere-tike 
2/3 < a/3 Disc-like 
2/3 > 2/3 Rod-like 
2/3 < 2/3 Blade-like 


Nichols “'°®) found that Y.» cannot be considered equivalent to 
y. An approximation @ to the previously described circularity 
¢ was defined as the ratio of the diameter of a circle of equal 
area to the diameter of the smallest circumscribing circle. Wadell 
felt that @ would be a good approximation to ‘Wop for roughly 
equi-axial particles only, and Pye and Pye“) confirmed it to 
be a poor approximation for other particles, since it is impossible 
to get a three-dimensional evaluation of an object from a two- 
dimensional measurement. 

Zingg proposed a rough shape classification, using the 
three mutually- -perpendicular axes a, b and ¢ (as eae 
defined for the ellipsoid) and these are employed as in Table 2 
From the results cited previously with regular shapes, this may 
prove to be a valuable supplement to the Y characterization. 
Cailleux “"® defined a flatness factor F, given by: 


given by: 


os 


and a roundness ratio R 


2 (Radius of sharpest corner) 
R, = aia (24) 


Longest diameter 


Von Andel, Wiggers and Maarleveld"") feel that y, the Wadell 
roundness and the Zingg critera are as easy to measure and 
more relevant than F and R.. Plots of Wop as a function of a/b 
and b/c have been calculated by Krumbein“!*®) who assumed 
ellipsoid shapes. 

Recently Aschenbrenner ‘'??) obtained an approximation to 
(which he called the working sphericity ~,4) in terms of a, b and 
¢ by assuming the particles to have the shape of a tetrakaide- 
cahedron. It is defined as 0.955 times the ratio of the surface 
area of an orthotetrakaidecahedron having the same volume as 
a tetrakaidecahedron, used for the approximation to the surface 
area of the tetrakaidecahedron. He found wW, to be a better 
approximation to Y than was Yos, especially at low W values. 


Krumbein “!®, Powers“! and Rittenhouse “'® have each 
offered very approximate visual means of characterizing particles 
by publishing illustrations of solids having various roundness 
and sphericity values. 

The significance of the various shape parameters discussed 
cannot be properly assessed until systematic drag experiments 
are performed with the many varieties of irregular shape which 
occur. It may very well be that some of the parameters which 
are extremely difficult to measure may not exert a noticeable 
influence in the overall momentum transfer. 


In all the investigations reported with freely moving bodies, 
the particles assumed an orientation which was dictated by the 
aerodynamical processes occurring in its boundary layer and 
wake. In actual multiparticle systems, frequent collisions with 
adjacent particles would Sommege the establishment of pre 
ferred orientation, and it would most likely become random. 
The same may be asi to result from free-stream velocity 
gradients or turbulence, this being supported by high- speed 
motion studies by Pasternak :'"'. Although it would require 
subtle experimental techniques, it is vitally important that drag 
determinations be made with single shapes under conditions 
which are more closely associated with the fluid parameters 
which are characteristic of multiparticle systems. ‘This is pres 
ently being attempted in these laboratories by Marchildon “°°, 
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Nomenclature 


Any consistent set of units may be employed. Those listed 

are merely illustrative. 

a = particle length perpendicular to largest circle or square 

cross-section, ft. 

shortest, intermediate and longest axis respectively of a 

triaxial particle, ft. 

A = projected area of sphere of volume equal to that of the 

particle, ft.? 

normal projected area relative to motion, ft.? 

area projected by a sphere with same surface area as 

particle, ft.? 

tangential area of particle relative to motion, ft.? 

Cp = total particle drag coefficient, dimensionless. 

total particle drag coefficient, dimensionless, Equation 

(10). 

Cp’ = empirical rotation drag coefficient, dimensionless, Equa- 

tion (1). 

drag coefficient for zero rotation, dimensionless. 

coefficient of lift, dimensionless, Equation (4). 

Becker “inertial drag’’ coefficient, dimensionless, Equa- 

tion (18). 

Cy = Becker ‘viscous drag"’ coefficient, dimensionless, Equa- 
tion (16). 

d = cylinder or sphere diameter, ft. 

d,, = arithmetic mean diameter, ft. 

d, = diameter of sphere of equivalent surface area, ft. 

d, = diameter of circle of same area as projection of particle 

d 

d 

d 


a,b,c 


ll ll 


S 
S 
& 

| 


e.¢ 
a ~] 
S 
til 


Qs 
ll 


lying in its most stable position, ft. 
min = minimum body length, ft. 
max = Maximum body length, ft. 
x = diameter of circle of area equal to projected area of 
particle, ft. 
diameter of sphere of equivalent volume, ft. 
= flatness factor = (6 + c)/(2a), dimensionless. 
Heywood spape factor = (2 /6)d,3/d.5, dimensionless. 
U,/U,;, dimensionless. 
Presler empirical velocity correction factor, dimensionless. 
K, = U,/U,,, for Stokesian region only, dimensionless. 
L = wetted diameter of particle, ft. 
mean hydraulic radius, ft. 
analytical coefficient, dimensionless. 
2 = analytical coefficient, dimensionless. 
sphere radius, ft. 
radius of edge roundness, ft. 
total drag force, Ib. force. 


2 (Radius of sharpest corner) 


a 
— 
Il 


Py 
| 


ll 


S55 
uuu 


ll 


> RLY: 
ll 


: = roundness ratio = - 
longest diameter 
dimensionless. 
Re = rg oh Number, dimensionless 
Re, = 4mp, U »/M, dimensionless. 
R, = Becker “inertial drag’’, dimensionless 
Ry = Becker ‘viscous drag", dimensionless. 
s = diameter of maximum cross-section or side of base of a 
square prism, ft. 
U, = mean velocity of particle relative to fluid, ft./sec 
U = velocity of sphere of equivalent volume calculated from 
Stokes’ Law, ft./sec. 
U, = peripheral velocity of particle surface relative to fluid, 
ft./sec. 
= volume of irregular particle, ft 
= downstream distance, ft 
= angle of incidence of flow, degrees 
= “normal factor” = A,/A, dimensionless 
w*r?/ L’,, dimensionless 
= tangential factor, Equation (11), dimensionless 
= time from beginning of motion, seconds 
= fluid viscosity, lb./(sec.) (ft.) 
= fluid density, Ib. /ft.4 
= particle density, lb./ft.® 
= longitudinal component of skin friction, lb. torce/tt.? 
= transverse component of skin friction, lb. torce/ft.? 
= form sphericity = md,?/4A,, dimensionless 
= circularity = (circumference of circle having the same 
cross-sectional area as irregular particle)/( perimeter of 
particle cross-section ). 
w = angular velocity, radians/sec 
tH = functional operator. 


3 


= 


PWDR R =~ - 
ll 


RENAE 


The Canadian Journal of Chemical Engineering, October, 1960 151 








U6/r, dimensionless. 

sphericity = (surface area of sphere of equal volume)/ 
(surface area of actual particle) 

(0.955) (surface area of an orthote trakaidecahedron 
having the same volume as a tetrakaidecahedron used 
for approximation )/(surface area of tetrakaidecahedron), 
dimensionless. 


Wil 


ll 


Yo.» = [(V)/(volume of smallest circumscribing sphere)]!'*. 
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Graphical-Numerical Determination of 


Variable Diffusion Coefhicient in Solids 


RAYMOND C. HALL? 


Considerable gain in information in characteris- 
tics of diffusion coefficient may be obtained if data 
are analysed for variable diffusion coefficient instead 
of for a constant value of diffusivity. This paper 
presents a simple method for graphical analysis of 
concentration gradients in solids to determine the 
diffusion coefficient as a function of concentration. 


Consideration is given to problems encountered in 
extrapolation of data. It is shown that, if relative 
values of diffusion coefficient is the objects of con- 
sideration, no great need exists for exacting extra- 
polation. Two extrapolations are employed on any 
given functional model; the results of one extrapola- 
tion are greater, while those of the other are less 
than the results predicted by the model. 


I" many papers dealing with diffusion coefficient the calculation 
of this coefficient involves at Geet an experimental determina- 
tion of a concentration gradient.“:?-5) “This gradient, usually 
shown as a graphical curve, is then “ie zed for the proportion- 
ality constant in Ficks’ law. In most real cases, especially as 
related to movements through solids, this proportionality factor 
or coefficient of diffusion is not a constant; frequently it varies 
with concentration. It is, however, common practice to assume 
this coefficient as a constant and to evaluate it, from the con- 
centration gradient, on this basis. 


Considerable information may be lost by such practice and 
one may wonder why data are not more commonly analyzed 
for variable coefficient of diffusion. It may be that many 
investigators have not followed through on methods of treatment 
of data in order to analyze for variable diffusivity. Crank“ 
and Jost“ present the basic equations by which a concentration 
gradient may be analyzed to coleuine D= f(c). Frequently 
the algebraic equation of C = f(x) is rather complex and 
difficult to write. In such a case, graphical, numerical or 
mechanical methods for integrations and differentiations may 
be employed. 

The objects of this paper are threefold: (a) to demonstrate 
the use of such graphical, numerical, and mechanical methods 
on an actual concentration gradient under conditions where D 
varies appreciably with concentration, (b) to point out and 
discuss several of the practical problems which were encountered 
in this solution, and (c) to demonstrate how these problems 
were handled. 

Although a specific solid, neat portland cement, and a specific 
moving material, solium, were used for this work, the methods 
employed can be applied to any suitable curve of concentration 
as a function of depth. Although the concentration gradient 
used in this work was established by means of radioactive 
tracer techniques, any method by which a suitable gradient can 
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be determined may be employed. Recent developments in 
techniques of use with radioactive materials have made these 
methods efficient and practical for experimental determination 
of concentration gradients in solids; hence this paper deals with 
a curve which was so established. 

In general, experimental procedures for studies on movement 
of materials in cement or concrete consist of preparing a test 
specimen having dimensions such that it can, for mathematical 
treatment, be considered as an infinite or semi-infinite volume. 
Actual dimensions of cement and similiar materials to meet this 
condition need not, in general, exceed about 1” square by 2” 
long with all the 1” X 2” plane faces sealed with an impermea- 
able coating. After sizing and sealing, the specimen is submerged 
in a prepared solution. The solvent of this solution is the 
material whose diffusion characteristics are under investigation 
and to which an appropriate activity of radioactive solvent has 
been added. After a suitable holding period the block is removed 
and activities at given depths are determined. From these data 
the geometric curve of concentration — depth is established; it 
is then analyzed for the diffusion coefficient. 

Prager“) developed equations by which the diffusion coefh- 
cient can be approximated as a function of concentration. For 
this method the total quantity of material moving into a semi- 
infinite film as a function of time is observed for a series of 
different driving potentials. Data are analyzed with a step by 
step method which requires a trial and error method on successive 
steps. This procedure has the advantage of being a non-destruc- 
tive method of obtaining data but does require data at several 
solution concentrations. The method presented below has the 
advantage of employing much simpler mathematics as well as 
requiring data at only one concentration of the solution phase. 

Recent advances in high speed mechanical computers and in 
programming for these computers have made the application of 
these devices to problems of this type very practical. As an 
initial step in writing a program for a digital computer, an 
analysis of a problem in a manner similar to the treatment 
given in this paper is frequently desirable. Hence this paper 
could serve as a basis upon which to write such a program. 

If only a few curves are to be analyzed — such as might be 
the case for a preliminary or for a feasibility study — it would 
be impractical to write a program. In this case the method 
outlined herein could be applied directly. 


Preparation of Concentration—Depth Curve 


It is assumed that the data of concentration as a function of 
depth is established. Normal variability in experimental data, 
inherent to all physical measurements, may be smoothed out by 
one of several methods. 


One method — readily adaptable to work in cement — was 
to plot the concentration on a logarithmic scale against depth 
on a linear scale. The best curve “by eye” was then drawn or, 
if suitable, a least squares line as described by Wylie“, could 
have been drawn. (Other approaches are suggested by Davis‘ 
and Snedecor"®).) The data of this best line were then re-plotted 
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Figui la. 


Evaluation of Diffusion Coefficient As a Function of 
Concentration 


Jost and Crank‘® suggest comparable methods for treat- 
ment of experimental data in order to evaluate a variable diffusion 
coeficient where the coefficient varies with concentration. The 
method involves integration and differentiation of the experi- 
mentally established gradient of concentration — depth. If the 
algebraic equation of this curve is difficult to write, methods of 
graphical or mechanical integration and differentiation may be 
empl ved. 

Using the equation, from either Crank or Jost, 


Cc. 
1 dx 
ee A eee dae (1) 


(, is generally assumed to be homogeneous throughout the 
block and usually is zero. 

It is assumed that concentration in the solution phase remains 
substantially constant and that all material which enters the 
block (the solid phase) does so by leaving the solution. This 
is comparable to providing a large quantity of dissolved solvent 
in solution relative to the total quantity of solvent which enters 
the test block. Equation (1) is comparable to 


AX oe Re 
), = — > MING? Bevescace sc 2 
E i ate 4) | rs c' Z\ ] ( ) 


where D; is evaluated at concentration C; on an actual concentra- 
tion gradient curve. The evaluation of (A X/AC); is shown in 
Figure la. 

A highly satisfactory means of establishing the perpendicular 
to the curve at depth x; is to employ the mirror method as 
described by Dana and Hillyard“. A front surfaced mirror, 
as provided by a polished metal surface, rather than a rear- 
surfaced glass mirror is required for this application. 


To draw the perpendicular to the curve at depth x;, the 
mirror is placed on the pivot C;, x; in such a manner that one face 
of the mirror is on a line which passes, perpendicular to the 
surface of the paper, through the pivot point and with the 
greatest length of the mirror approximately perpendicular to 
the curve. The mirror is then rotated on this vertical line 
which passes through the pivot point until the visual reflection 
of the curve and the curve on the paper join so as to form a 
smooth and continuous curve. The position of the ends of the 
mirror are then marked on the paper. The procedure is next 
repeated, about this same point, on the other segment of the 
curve. The perpendicular to the curve is then drawn through 
C;, x; such that it passes (as best possible) through the points 
where the ends of the mirror were marked and the perpendicular 
to this line is then drawn. The slope of the last line which is 
parallel to the tangent is taken, after proper scaling as AC/AX. 
Or the points which marked the ends of the mirror may be 
employed directly, with attention to signs and scaling, to make 
this calculation. 


One method for evaluating 


> XAc 


is illustrated in Figure 1b. 


In this case the trapezoid rule is followed. 


C=C, Xo+X ¥.4+X 
zy Xx t-t.i—— G=-G 1 
Bg wace(ona) (™F*)+(4-4)(*5")s 
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Figure la—Use of mechanical and graphical methods to 

evaluate AX/AC at X,. A front surfaced mirror is employed 

to establish the perpendicular to the concentration gradient 
curve at point X,, C,. 

Figure 1b—Use of the trapezoidal rule to evaluate Equation 


(3). Region of questionable extrapolation, under dotted 
curve = (C, — C,)(X, — X,)(%). 


If greater accuracy is required, methods such as Simpson’s 
rule“), or improvements by extrapolation as outlined by 
Salvadori and Baron“), may be employed. 


The evaluation of 1/2t is straightforward. 


Thus, by repeating this procedure at various X, and calculat- 
ing the corresponding D;, the plot of D; against C may be 
prepared. 


Sources of Difficulty 


It is observed that the first term in the graphical integration 
by the trapezoidal rule necessitates an evaluation of X, where 
the curve conv erges to C,. Various algebraic expressions which 
may serve to fit experimental data in a satisfactory manner, 
such as equilateral hyperbola or reciprocal exponential type, 
may not converge to C = 0 for X<@. This, however, presents 
no real problem i in obtaining suitable accuracy in evaluating 


€ = G 
2 XAC 
C=C, 


The effective value for _ may be set at any value for C sufh- 
ciently small and the AC = (C, — C,) taken to any practical 
degree of smallness. Such procedure 3 is justifiable because, by 
the nature of the system and the treatment to which it may 
have been subjected for any finite t, the actual curve must and 
does converge to C, at finite x. Failure of an algebraic equation 
to converge to zero is merely an example of how an idealized 
algebraic equation can fail to adequately describe the real 
physical situation. Thus, two major problems develop; the 
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first is the evaluation of the concentration gradient at depths 
beyond which measurements were made and the second is the 
evaluation of the initial concentration of the moving material in 


the solid. 


In order to demonstrate these, as well as other practical 
problems which were encountered in ev aluating the integral and 
the differential at various points, some data of Chang 4) were 
analyzed. The techniques applied are presented here. 


The Problem of Extrapolation 


The depth — concentration data were assembled in Table 1*. 
The block of concrete (2” diameter cylinder, +” long, and having 
all surfaces except the ends covered with wax) was held sub- 
merged in a radioactive solution for a period of 273 days. 

The two ends of the block, Ends I and II, were analyzed 
192 days and 635 days respectively after removal of the block 
from the radioactive solution. It was reasonable to assume that 
both ends of the block were of similar physical and chemical 
characteristics, so that diffusion characteristics should have 
been very similar. It is obvious that the concentration depth 
curves for the two ends do not coincide. Although it is not the 
purpose of this paper to discuss the reasons for this discrepancy, 
it may be noted that the total sodium in the block (i.e. the 
integrated area under the curve) was greater for the measure- 
ments which were made on End II than for those made on End I. 
Since the block was not exposed to any external source of 
sodium during this time it appears as if av iol: ation of the material 
balance was obtained. Since this may be taken as an error the 
question was naturally raised as to w vhich curve was represent- 
ative of the true situation. One suggestion was that both curves, 
used together, were the most representative of the true data, 

Data were plotted on semi-log paper as shown in Figure 2 
where plots for counts, as made on Ends I and II, 92 days and 
635 days respectively after removal of the block from the solu- 
tion, are plotted as separate curves and the least squares curve 
for both sets of data combined together is also shown as one 
straight line. It was noted that the best curves by eye tended 
toward an § shape which was distinctly different from a straight 
line. Data from each of these three curves were plotted on 
rectangular coordinates. The plot from the resulting “best by 
eye’ curve for the data of End II is shown 1n Figure 3; it Is 
generally representative of the other two curves. 

Steps 1, 2, 
the mirror method to the evaluation of 4, X//A\C. 

The problem of extrapolation of the curve in Figure 3 was 
considered; this extension was necessary in order to effect 
the integration of (X/.C) from C, to C;. 
the curve on Figure 2 was made. Extrapolation of the least 
squares line on this plot was easy; 
how closely this line came to true values was raised. 
decided to graphically extrapolate the “best by eye” 


It was 
curve. 

‘The practice of using extrapolations is usually considered as 
hazardous. Even if a suitable plot or mathematical model can 
be found which will give “‘best fit” over reasonable ranges of 
data there is no guarantee that this straight line or this set of 
evaluated constants will remain valid over appreciable extrapola- 
tions. Obviously it is desirable to obtain data on concentrations 
at appreciable depths in order that the severity of this source of 
errors may be reduced to a minimum, Accurate evaluation of 
concentrations at low values is often difficult and frequently is 
costly to make. Hence, some method of extrapolation will 
normally be required 

\ common type of study involves consideration of a large 
number of specimens all having essentially the same fundamental 


*Tables 1, 2, 3, and 4a of this paper have been deposited as Document 

No. 6394 with the ADI Auxiliary Publications Project, Photoduplication 
Service, Library of Congress, Washington 25, D.C. A copy may be ob- 
tained by citing the Document No. and by remitting $1.25 for photo- 
prints, or $1.25 for 35 mm. microfilm. Advance payment is required. 
Make cheques or money orders payable to: Chief, Photoduplication 
Service, Library of Congress. 
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and 3 on Figure 3 illustrate the application of 


‘The extrapolation of 


however, the question of 


characteristics. As an example, one may be concerned with the 
influence of density of material on the D values. A series of 
specimens would be prepared, all having different densities. If 
suitable evidence is available, or may be established, to indicate 
that the general characteristics of the concentration gradients 
are independent of density (even though slopes or constants of 
a mathematical model may vary between densities), one would 
be justified in expending great effort and expense to obtain data 
at low concentrations on one speciman in order to establish the 
general mathematical model. Data on the remaining specimens 
would then be obtained at only enough points to evaluate the 
constants of the model and these would then be used to extrapo- 
late to low C. 

In other cases sufficient information, either experimental or 
based on postulated mechanisms, may be available to aid in 
writing the mathematical model. 

In the event that no information is available and expense or 
other restrictions prohibit its establishment, an extrapolation 
will be necessary. This paper suggests that it is not necessary 
to obtain an exacting extrapolation: rather it proposes that the 
two extrapolations be made; one having values greater than and 
one having values less than the apparent true value. The 


C=n 
resulting calculations will provide two sets of values (of 2 xde), 
C=0 


one of which may be assured to be less than true and one greater 
than true. The range of these values of D may then be considered. 
If the range is too great for the desired use, an improved means 
of extrapolation must be found. If, on the other hand, the range 
is within reasonable limits for the immediate purpose at hand, 
the extrapolation method used may be considered as satisfactory. 


The curve of Figure 2 is one which must be extrapolated. 
The data were taken several years ago by another investigator 
and no method was available by w hich data at greater depths 
could be obtained. The specimens were several years old at 
the time they were submerged and they remained in solution 
for nearly a year; it would take a long treatment time to repro- 
duce the work. No information was available on which to base 
a hypothetical mathematical model. 


A first method of extrapolation is demonstrated in Figure 4 
Data had initially been smoothed on a plot of log C = f(x), 
(Figure 2). The end of this curve, on the semi-log paper, was 
plotted on Figure 4. On this plot it appeared as if the end of 
the concentration gradient was approaching a straight line. Two 
extrapolations, both straight lines, were made. E xtrapolation A 
was made with slope appreciably greater and extrapolation B 
was made with a slope appreciably less than the apparent slope 
of the end of the data curve. 

A second approach to extrapolation is demonstrated in 
Figure 5. In this case the (C/A X was plotted against depth. 
The resulting plot indicates that the slope is approaching a 
constant value at the greater depths. From this curve two 
values of slope were chosen for extrapolation, one greater and 
one less than the apparent final slope. 

In a third method, demonstrated in Figure 6, a model of 
the form C = A*~ was assumed. For depths in excess of 0.07 
inches it appeared as if the values of A and b were constant. 
The resulting straight line was extended to depths of 30 inches. 
Only one line was employ ed in this case although two, one a 
maximum and one a minimum, could have been used. 


Calculations and results of the extrapolations A and B from 
Figure 4 a in conpenetinn with the curves from Figure 2 are 
shown in Tables 2 and 3* and are plotted in Figure 7. Compar- 
isons of the values of D as obtained from the various extrapola- 
tions of Figures 4, 5 and 6 are tabulated Table 4a* and 
plotted in Figure 8. 

In as much as the extrapolation of Figure 6 extends the C 
to very low values and, at higher x values, causes the curve to 
be nearly parallel to the x axis, the need for adequate evaluation 
of C, was brought into focus. This is demonstrated by con- 
sidering three different values for C,. 
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Figure 2—Concentration gradient of sodium in concrete 
specimen. 
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Figure 4—Extrapolations (for End II of curve by eye from 

Figure 2) to C = C,. Extrapolation B at slope less than 

apparent slope of data curve and extrapolation A at slope 
greater than apparent slope of the curve from data. 


The Problem of Initial Concentration 


The importance of reasonable evaluation of C, is involved in 


Cy 
the evaluation of  X/\C where 
Os 
Cc. Ge Ge 
ZXAC = TXAC+ SXAC. (4) 
C; C; C, 
when x, is the deepest point where /(\X//\C is measured and 
x; is the point where Dj is being evaluated. Thus if 
Ge 
wo ARC = Ee... (5) 
Ge 
then 
ic. Cc. 
DAAC 2 XAC+K (6) 
c. G 


and the problem is essentially that of evaluating K where the 
magnitude of K depends on the value of C,. In actual practice 
many natural materials may contain a me: somrediile degree of 
radioactivity which may be of importance in the evaluation of 
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Figure 3—Concentration gradient for End II. “Best” curve 
by eye from Figure 2 replotted on rectangular coordinates 
and analysed for (AX/AC,). Extrapolations with constant 
slopes of 300 and 800 (Figure 5) also shown in dotted 


lines. 
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Figure 5—Plot of slope at various depths (for End II of 

curve by eye from Figure 2). Assume extrapolated curve 

holds constant slope of between 800 and 300 counts per 
minute/inch. 


C; 
~ NAC, especially if introduced activities are relatively low 
C 


or if extrapolations, as mentioned above, are critical. In this 
particular study the evaluation of C, was inadvertently neglected. 
The potential seriousness of this error was further considered 
by assuming a natural background of C, = 20 ¢/m on data 
obtained 92 days after removal of the block from solution. 
These calculated values were compared with those which were 
determined on the basis of no natural background. Results are 
summarized in Figure 9 (extrapolation of the curve by eye was 
made by assuming a slope the same as for extrapolation B of 
the data obtained 635 days after removal of the block from the 
solution). 


Conclusions 


Results of the application of the proposed method of analyz- 
ing for variable diffusion coefhcient to the system sodium- 
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Figure 7—Comparison of D values based on extrapolation 
of least squares line and extrapolations A and B (Figure 4) 
of data from Figure 2. 


cement, clearly demonstrated that the coefficient is not a con- 
stant but varies appreciably with concentration, especially over 
the lower ranges. 

Results of various methods of extrapolation demonstrated 
that, if information on only relative values of D be needed, no 
great concern need be given to the method of extrapolation 
(Table 4b). If, however, exacting values of D are required, a 
method must be found which will predict the shape of the 
concentration gradient to considerable depth with considerable 
accuracy. If no other means are available upon a graphical or 
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Figure 8—Comparison of D values from Table 4a. Extrapo- 

lations A and B from Figure 4. Extrapolations 800 and 300 

from max. and min. of Figure 5. Extrapolations C, = 20, 
5 and 2.83 from Figure 6. 


a mathematical model can be built to accomplish this extrapola- 
tion, it is suggested: that data necessary to construct a complete 
curve be gathered on at least one representative specimen, that 
the mathematical model be built from this information and that 
this model be used for extrapolation of fewer data on brother 
specimens. 

In the problem considered there was no alternative to making 
an extrapolation. The model, C = At gave D values having 
order of magnitude comparable to D values which were calcu- 
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lated by assuming a constant coefficient. It was felt that this 
method was superior to others, if true values of D were desired. 


Need for evaluation of inherent concentration was demon- 
strated, if true value of D is needed. 


Wider use of a method such as this is suggested as a means 
of obtaining more and better information on D values by analysis 
of concentration gradients. 


Nomenclature 

C = concentration, convenient units. 

C, = initial concentration in block 

C; = concentration in block after time ¢ and at depth Xj. 

D = coefficient of diffusion, length squared per time 

D; = coefficient of diffusion at a given C;. 

t = time interval during which block was submerged. 

XY = length (depth into block, measured from and perpendic- 
ular to exposed surface) 
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Concentration Distribution in 


Two-Phase Pipe Flow’ 


S. L. SOO? and J. A. REGALBUTO* 


Concentration distribution of solid particles in 
two-phase turbulent pipe flow was studied. Due to 
the inertia of solid particles, transport of solid par- 
ticles across the velocity field of the turbulent stream 
by turbulent diffusion results in deceleration of par- 
ticles in the region near the wall by the stream and 
acceleration of particles transported toward the high 
velocity core. In the core, the particles have lower 
mean velocity than the stream while near the wall 
the mean particle velocity is higher. This gives rise 
to a non-uniform distribution of solid particles in 
fully developed turbulent pipe flow. Experimental 
determination of concentration distribution was made 
with an impact disk particle counter system. 


Ee tee of wall interference in two-phase (solid particle- 


yas system) turbulent motion indicates that the presence of 


the wall such as in fully develo »ped pipe flow affects considerably 
the mean stream condition’:?). The presence of the wall in- 
creases significantly the intensity of particle motion above that 
in free turbulence due to continuous transport of momentum 
from the wall to the core. It is reasonable to question the effect 
of momentum transport from the wall on the concentration 
distribution over the duct. 


\ knowledge of the distribution of the solid particles 1s 
significant, considering its effect on heat and mass transfer 
processes usually encountered in the applications of two-phase 
turbulent motion. In References and * uniform distribution 
of the concentration of solid particles was assumed. The experi- 
mental results of Reference considers basically the gravity 


effect on two- phase streams of large weight ratio (loading) of 


solid particles to gas, and when there is no apparent gravity 
effect, the concentraticn of solid particles was taken to be 
uniform. This paper concerns the concentration distribution in 
the case of smal] loading, where gravity effect is measured to 
be small. 

At present, very little is known about pipe turbulence of a 
single phase fluid. Various conjectures were made as a matter 
of qualitative explanation of experimental results of mean veloc- 
ity distribution, turbulence intensities, shear stresses, eddy 
viscosities, turbulence dissipation, energy balances, and energy 
spectra. A plausible explanation ts that there is a flow of kinetic 
energy by turbulent diffusion down the gradient of turbulent 
energ toward the cenver of the pipe where eddies are dissipated, 


Manuscript received December 27, 1959; accepted July 15, 1960 
Professor of Mechanical Engineering, University of Illinois, Urbana, Ill 
Texti'e Research Fellow, Princeton University, Princeton, N.J. His con- 
tribution *a the experimental program was in partial fulfillment of 
requirements of M.S. degree 

Based on a paper presented to the C.1.C. Chemical Engineering Con- 
ference, November 9-11, 1959 


160 


while turbulence diffusion of pressure energy is directed from 
the core region toward the wall. There is continuous transport 
and diffusion of small energy-rich eddies from the wall region 
toward the center. At the same time there is a flow of energy- 
poor fluid lumps back toward the wall, where mean motion 
shears these lumps into small high-intensity eddies. The 
phenomenon is extremely complicated and comprehensive ana- 
lytical treatment is rendered difficult. The net situation is 
that there is turbulent transport of momentum through accel- 
erating the fluid from the wall toward the core region. 

It might even be said that our knowledge of the one-phase 
pipe turbulence has not yet readied us for studying a two-phase 
solid-gas suspension. The significance of the problem itself, 
however, calls for some immediate attention. Due to transport 
of fluid momentum toward the core region, as explained in the 
above, the solid particles are transported by the fluid, opposing 
the diffusion mechanism due to concentration (number density 
of solid particles) gradient. It can be expected that the distri- 
bution of solid particles i in fully developed two-phase turbulent 
motion will not be uniform over the cross-section of a pipe. 


In addition, the inertia of the solid particles play an important 
part in the transport process. Solid particles transported from 
the center to near the wall is decelerated to the local mean 
stream velocity. The mean velocity of solid particles uy» in 
this region 1s therefore higher than that of the stream uw. Based 
on the local static number density or concentration, C, (number 
of particles per unit volume of fluid when both are at rest or 
both moving at the same speed), the local number density C 
is then 


Riis es (1) 


\ higher particle velocity than the stream velocity shows up as 
a decrease in particle concentration. 

On the other hand, the solid particles transported from near 
the wall (which are already decelerated) to the core region is 
accelerated by the stream to the mean velocity of the core 
region. The average particle — u, in the core region is 
lower than that of the fluid, Based on Equation (1), this 
shows up as an increase in ania concentration above the local 
static value. 

The net effect of transport of momentum and inertia of 
solid particles (especially the latter) can be visualized, without 
going into the details of particle motion itself, as acceleration 
of solid particles transported by turbulent diffusion from the 
wall toward the center of the pipe and deceleration of solid 
particles transported from the center toward the wall. The 
direction of acceleration is therefore radially inward toward 
raising the particle concentration in the core, opposing the diffu- 
sion mechanism due to concentration gradient. It should be 
noted that one cannot account for the effect of particle inertia 
by considering any special distribution of turbulent diffusivity 
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Figure 1—Schematic diagram of apparatus. 


of particles alone; it would give uniform distribution of solid 
particles when fully dev eloped turbulence is achieved. 

The above phenomenon can be formulated in the following 
manner as a first approximation. 

The rate of radial motion of solid particles due to concentra- 
tion gradient is given by: 


n= —D (2) 
dr 


where r is the radial coordinate, D the particle diffusivity (due 
to stream diffusivity)’, C the particle concentration. 

Neglecting gravity effect, the change of particle concentra- 
tion by diffusion, or the total variation of u with respect to 
time is: 


dn a u Oo D le ; 
= (rn) = r (9) 
dt r Or r OF 0; 
for steady motion (0/0t = 0) at mean axial flow velocity uw. 


When a steady concentration profile is maintained, the rate 
dn/dt is due to transport of — across the fluid velocity 
field and is proportional to the drag force acting on the particles. 
Hence as an approximation, 47 dt is proportic nal to the concen- 
tration C at radius r, or 

dn 


= fC, (4) 
dt 


f is the acceleration of solid particles transported by turbulent 
diffusion from the wall toward the center of the pipe and decel- 
eration of solid particles transported from the center toward 
the wall. f relates fluid drag and inertia of solid particles as a 
solid particle moves across the distributed vel city field of the 
turbulent stream. 

Combining Equations (3) and (4), we get, for constant 
particle diffusivity, 


DO ac fC 4 
ao + = () (3) 
r Or 0) Mu 


f, in general, is a function of state parameters of the system 


and radius r. 


If an average value of f is taken, the solution of Equation 
(5) is then: 
C= C,J, (ry f/Du), (6) 


where C, is the concentration at the center of the pipe or the 
reference concentration, and /, 1s the zero order Bessel’s function. 
The above relation gives a concentration profile somewhat simi- 
lar to that of mean velocitv. The parameter affecting concen- 
tration distribution /R? Du (where R is the pipe radius) can be 
determined from measured concentration distribution. 

\n impact disk particle counter was developed for the pur- 
pose of measuring the distribution of solid particles in two-phase 
turbulent pipe flow (glass beads in air). Measurements were 
made on a two-phase flow facility. Figure | shows the duct, 
the solid particle sampling system and the particle counting 
system schematically. 


Two Phase Flow Facility 


The measurements of concentration distribution were made 
along a vertical diameter of a three-inch inside diameter pipe. 
Air is induced through the pipe by a Buffalo Vol:me Fan driven 
by a 5 hp. shunt- wound D.C. motor. ¢ ‘hanges in the air stream 
velocity are facilitated by field control of the motor or manipula- 
tion of a bypass valve immediately preceding the blower intake. 
The dimension as shown in Figure | ts such that fully developed 
pipe flow exists at the points of sam ipling”. 

Foreign particles are excluded from the air stream by approxi- 
mately 15 sq. ft. of § micron filter paper surroun ling the pipe 
intake. A pitot-static tube located 10 diameters downstream is 
used to report the mean stream velocity via an inclined draft gage. 

The solid phase (spherical glass beads) 1s metered into the 
air stream 20 diameters downstream from the air inlet by an 
inclined piston and cylinder arrangement. The uppermost end 
of the cylinder is connected to the main pipe by a short vertical 
tube. Beads loaded into the evlinder are torced upward by the 
piston over a shallow V-notch weir, thence down the vertical 
tube into the air stream. Linear motion of the piston is accom- 
plished by a power screw and nut mechanism. The nut is 
secured to the piston and constrained from rotating while the 
screw is powered by a | 3 hp. D.C. motor driving through a 
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Figure 2—Schematic diagram of sampling system. 


Zero-Max variable speed transmission. This combination pro- 
vides both the torque necessary to rotate the screw and the 


means by which to regulate the screw r.p.m., and hence the 
bead feed rate. 


A theoretical feed rate, calculated trom power screw r.p.m., 
thread pitch and the cylinder diameter compared favorably with 
that obtained by w cighing a sample collected by the separator 
over a measured time interval. Further calibration was deemed 
unnecessary and the tachometer indicating the screw r.p.m. was 
calibrated to give the feed rate directly in Ib./hr. 

At the downstream end of the pipe, the glass beads are 
divorced from the air by diffusing the stream into a cyclone- 
type separator. The air exits through 3 sq. ft. of 5 micron 
paper at the top of the separator w hile at the bottom, a double 
door device allows collection of the beads without unduly 
disturbing the equilibrium of the flow system. 


Sampling of the Solid Phase 


The solid phase distribution is sampled by a probe fabricated 
trom .062” inside diameter, .016” wall tubing. The probe is 
mounted on a vernier depth gage which allows easy traverse in 
the duct and further, locates the probe precisely (to within .001 


inch). 
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A 90°, 1/2 inch radius bend in the tube collects the bead 
sample downward out of the main duct into a vacuum bell 
housing the transducer of the counting system. A collecting 
cone then funnels the sample into either of a pair of weighing 
bottles, after the beads have collided with the transducer. 
(Figure 2) 

To insure that the sampling is truly representative of the 
actual particle distribution, the vacuum bell is kept at a pressure 
low enough to induce, through the probe, an airflow whose 
velocity at the probe inlet is the same as the air stream velocity 
in the immediate vicinity. In the line between the vacuum 
pump and the bell, a bypass valve (opening to the atmosphere) 
controls the air flow through the probe. Thus the velocity 
within the probe can be matched not only to different average 
air velocities but, more important, to the velocity profile en- 
countered as the probe traverses the duct. 


A static pressure tap, located 8 diameters back from the 
probe inlet in a horizontal plane, is connected by a 1/32 inch 
inside diameter tube to one leg of a U-tube manometer. The 
other leg of the manometer is connected to a second static tap 
in the wall of the main duct at the same station as the probe 
inlet. Since, for the least disturbance to the flow field, the veloc- 
ities within and without the probe are desired to be equal, it is 
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necessary that the static pressures measured by the manometer 
be identical. Thus the velocities can be matched simply by 
adjusting the vacuum bell bypass valve until the legs of the 
U-tube manometer are equal. An added feature of the mano- 
meter is its adjustable inclination which allows matching of the 
velocities to any resonable degree of accuracy. 

Finally, added protection is given the vacuum pump against 
ingesting glass beads by the design of the aspirator shown in 
Figure 2. Particle inertia is utilized to separate the solid and 
gas phases as the air stream must reverse direction 180° and 
then pass through a closely packed steel wool filter before 
reaching the vacuum pump. 

Substantiation of the velocity matching capability of the 
induction was first attempted by photographing the two-phase 
flow in the vicinity of the probe. The photographic technique 
and equipment used in this test is described in Reference”). A 
series of pictures was taken at various values of the difference 
in static pressure (A\p,) within and without the probe. Checks 
were made with flow conditions at 40 ft./sec., 115 micron 
beads (more susceptible to suction of probe than 230 micron 
particles) and average loading of .025 lb. of beads per Ib. of 
air. Photographs* show that the bead flow pattern in general 
was not changed significantly by the presence of the “probe, 
The gas flow pattern might be affected locally, but we are only 
concerned with motion of the beads. This is shown by the near 
straight line trajectories in all the photographs referred to above. 

However, the photographs could not give enough detail 
about the effect of Ap, on the quantity of beads ingested by 
the probe. To this end a second test was made which resulted 
in a plot of beads ingested (in terms of the count recorded per 
second) versus A\p, (Figure 3). 

The shallow slope at Ap, = 0, indicates that the number of 
beads ingested is relatively constant if the manometer legs are 
equalized within 0.1 inch total. 


Actually no accurate adjustment for /\p, is necessary as 
long as it is constant for a given set of traverse. As shown in 
Figure 3, within the range where A\p, vs. counts per second 
is linear, when the number is directly proportional to Ap,, 
subsequent normalization rather than actual number of count 
determines the concentration distribution. This normalizing 
procedure is necessary also because variation in particle sizes 
within the range of 210 to 250 microns for 230-micron particle 
diameter and 105 to 125 microns for 115-micron particle 
diameter. 

To insure that samplings were made at points of fully devel- 
oped pipe flow, checks were made at locations 3 feet apart along 
the pipe. The concentration distribution was shown to be steady. 


Solid Phase Counting System 

An ordinary Sonotone 3T ceramic cartridge, intended pri- 
marily for use in high fidelity phonographs, is the transducer of 
the “counting” system. The cartridge 1s modified by replacing 
the needle with an elliptical plate connected to the ceramic 
crystal by a short stem. The axis of the stem is oriented in 
such a way that beads impacting on the plate excite the crystal 
in exactly the same manner as the original needle. (Figure 2) 


Anticipation of a low signal level from the crystal dictated 
the use of a low-noise, high-gain, battery-powered amplifier. 
Chosen was a Tektronix Model 122 with a maximum gain of 
1000 and a maximum noise level of 20 microvolts r.m.s. 


After amplification, the signal is fed into a Hewlett-Packard 
Model 521A electronic counter and an oscilloscope. A useful 
feature of the counter provides control of the input voltage level 
to the counting circuit. Thus, given a favorable signal-to-noise 





*Figures 3a, 3b, 3c, 5b, 5c, 5d, 6b, 6c, 7b, 7c, 7d, and 7e, of this paper 
have been deposited as Document No. 6395 with the ADI Auxiliary 
Publications Project, —— Service, Library of Congress, 
Washington 25, D.C. A copy may be obtained by citing the Document 
No. and by remitting $1.25 for photoprints, or $1.25 for 35 mm. micro- 
film. Advance payment is required. Make cheques or money orders 
payable to: Chief, Photoduplication Service, Library of Congress. 
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Figure 3—Counting system response to variation in | p,. 


ratio, unwanted noise and random vibrations can be excluded 
from the count. Precise time gates of one-tenth second and one 
second are provided along with a manual gate for counting over 
an extended period of time. 

General information about signal-to-noise ratio, signal wave 
form, presence of a signal, etc., is checked by an oscilloscope. 

Undesirable interference from mechanical or electrical 
sources is eliminated insofar as possible. The transducer is 
isolated from mechanical vibration by extensive foam rubber 
mounting. Connecting tubes and cables are as flexible as possible. 
To prevent stray electrical pickup, the cases of the various 
components were connected and the system grounded at a 
single point to prevent circulating currents. Another important 
item in this system is a wire mesh cage enveloping the transducer. 


Transducer Calibration 

The response of the modified cartridge was investigated by 
subjecting the impact plate to sinusoidal air pressure variations 
and view ing the output on an oscilloscope. To obtain the pres- 
sure variations, the output from a variable frequency oscillator, 
suitably amplified, was broadcast from a high fidelity speaker 
system. Varying the oscillator frequency showed resonances 
at 400 c.p.s. and 800 ¢ c.p.s. in addition to a response limit of 
approximately 8000 C.p.s. The resonance amplitudes were 
small enough in comparison to the bead signal and the frequency 
response high enough in comparison to the anticipated average 
number of counts per second, so that further modifications were 
deemed unnecessary. 

An initial calibration to allow conversion from counts per 
second to pounds of beads per second was made by weighing a 

“counted” sample collected under operating conditions. The 
count was manually started at the instant the two-way valve in 
the collecting system was switched from the stand-by bottle to 
the weighing bottle. 

Weighing inconsistencies, while not unduly large, resulted 
from the extremely small samples obtained over a long time 
period. Thus an alternate method utilizing the average bead 
flow was used. Plotting each solid phase distribution profile (in 
terms of counts per second), a mean effective figure of counts 
per second was obtained. Comparing this average figure with 
the known average bead flow indicates the necessary magnitude 
of the desired conversion factor in pounds per count. This 
latter method is preferred because of the larger quantities in- 
volved and the greater accuracy to which these quantities are 
known. 

Due to the lack of a precise time gate longer than one second 
in this particular counter, it became necessary to determine the 
number of intermittent digit displays to be av eraged for a repre- 
sentative bead flow figure in counts per second. 
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Figure 4—Velocity profiles at probe location. 


Averaging a 60 second sample after every 10 display s showed 
that the average of approximately 20 displays resulted in a figure 
only 1% different from the overall average. 


Determination of Concentration Distribution of the 


Solid Phase 


At each combination of average stream velocity, bead size 
and loading, bead flows were measured in one-tenth inch incre- 
ments from wall to wall. 
was allowed to equilibrate and final equilibrium” was considered 
attained only after maximum and consistent counts per second 
were observed. 


Prior to recording the 20 consecutive displays at each point 
in the traverse, the only adjustment necessary 
matching velocities within and without the probe. 
the feed rate and average stream velocity 
adjusted if necessary. 


However, 
were checked and 


In this manner fluctuations were minimized to approximately 
0.1% in the feed rate and 2% in average stream velocities. 


Results 


The original data used to obtain the concentration distribu- 
tion is itself of interest and is presented graphically. The 
velocity profiles, from which the air flow at a discrete point 
is determined, are given in Figure 4 at average air (at room 
temperature and atmospheric pressure) velocities of 50, 75 and 
100 feet per second. 

Solid particles of two mean diameters 230 microns and 115 
microns were used. The loading chosen were .05, .10, and .15 
lb. of solid per Ib. of air. 

Figure 5, 5(b), 5(c) and 5(d) give the counts per second 
(N) obtained by the solid particle counter. The mass flow is 
obtained by determining a normalizing constant K, in such a 
way that , 

(CR 
K 2x| Nrdr = 


o 


Total solid flow rate (7) 


and the mass flow of solid particles is obtained from the product 
NK, where K may be expressed in the unit lb./in.? per sec. per 


164 


Before recording any data the flow: 


was that of 


O 05 LB SOLD 
4 10 


LB. AIR Ay. 


oO 158 


75 FPS STREAM 
230 MICRON SOLIO 


IN 
° 


RADIUS, 





oO 50 100 150 200 250 


COUNTS /SEC 


Figure 5—Variation of counts with position. 


(c.p.s.), giving the mass flow in the unit of |b./in.? sec. or Ib./hr. 
in.2. Figure 6, 6(b) and 6(c) show results calculated for the 
case of 230 micron glass beads in air. 

The concentration of solid particles is obtained by dividing 
the local rate of mass flow of solids by the local mass flow rate 
of air based on Figure 4. Figure 7, 7 7(b) and 7 (c) give the con- 
centration distributions in the case of 230-micron glass beads. 
Figure 7(d) and 7(e) give the mass flow and concentration 
distribution of 115-micron glass beads in air. 

The parameter fR?/Du can be determined from these meas- 
ured data. Since, from Equation (6), 


R (rR y 
2x | Crdr = on, | J, a—)rdr. (8) 
oe o R 
where a = ¥fR?/Du 
if weletar/R = 
bate ie 
a | JAX) XAX = my Gleick... ae 


The right-hand side of the above equation can be calculated 
from experimental data. The integral on the left hand side is 
given in Figure 8, from which a or f{R°/Du can be determined 
by curve fitting as shown in Figure 7 

The calculated results are summarized in Table 1 


Discussion 

It is evident that an interesting regime in two phase turbulent 
How has been encountered by the combinations of particle size, 
average velocity, and loading chosen for this investigation. The 
concentration profiles clearly indicate that the effect of gravity 
on the particles becomes s negligible (1) as average stream velocity 
increases (average loading remaining constant) and (2) as the 
bead size decreases. However, the concentration profiles at 
100 f.p.s. indicate that increasing the loading at this higher 
velocity makes the gravity effect less significant. This suggests 
the existence of particle interaction) even though the mean 
distance between spheres at a loading of 0.15 Ib. of 230 micron 
beads per Ib. of air is approximately 20 diameters. 
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The most striking common feature of the concentration 
profiles is the maximum which occurs (in some cases, spectac- 
ularly) at the center of the duct. 

The investigation stopped at an average stream velocity of 
50 f.p.s. and 115 micron beads because beyond this point (higher 
velocities and smaller bead sizes) the signal-to-noise ratio rapidly 
approached one. That is, the pressure fluctuations of the turb- 
ulent air stream impinging on the transducer plate produced 
signals indistinguishable, except by form, from the bead signals. 


Any further investigation would necessitate the addition of 
a discriminating circuit between the transducer and the counter. x 
This circuit would be designed to pass the short rise time 
bead signals while rejecting the approximately sinusoidal noise 
signals?) 


Obviously, the value of f should decrease to zero toward = 
the center of the pipe and w increases toward the center of the ~ 
Pipe. Table 1 shows the average values of f/R?/Du and f of te. 
various cases. Rather than making detailed inferences regarding —% 


the trend of variation, it can be said that the measurements, 

within its range of accuracy, suggest a value of the order of 2 

for the wall effect parameter /R? Du in the case of 230 micron 6 P oO . 7 r. 
particles and 4 for 115 micron particles. This is for a pipe , : = 
flow Reynolds number range of 80,000 to 140,000 and particle a -/[# 

diffusivity, D, at the ratio of density of solid to density of gas 

of 1780, (density of glass/density of air). The acceleration f Figure 8—Diagrams for the determination of concentra- 
is naturally smaller for lower stream velocity. However, for tion distribution. 
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TABLE I 
Summary of Results 























Solid Mean i a Loading Solid 
particle stream ae ee Ib. solid | flow fRIY | f | f/u 
diameter, velocity e "te tVaee - per rate / Du |  f.p.s.? | 1/sec. 
micron | f.p.s. | ear Ib. air | lb./hr. 
a ——_ ze | 5 
230 100 2.7 x 10-* | SiS | 98.4 | 2.19 38.0 .380 
| .10 | 65.5 4.16 ten 721 
.05 | 32.8 | 4.36 75.7 .757 
75 2.4.x 10-8 15 | 147.3 3.45 39.7 .529 
10 | 98.2 3.38 39.0 .520 
.05 | 49.1 4.03 46.5 .620 
50 2.0 x 10-8 1S 196.5 3.67 23.5 470 
.10 131.0 4.00 25:6 «oz 
.05 65.5 1.69 10.8 .216 
115 50 1.4x 10" 55 98.4 4.16 18.0 . 360 
.10 65.5 4.94 21.4 .428 
.05 32.8 4.40 19.0 . 380 
different particle sizes involved, no significant difference was Nomenclature 
found for the same mean stream velocity. Therefore, the value C = Concentration or number density of solid particles (or 


of {R*/Du depends to a large extent on the value of D. The 
larger the value of D (very small particles, high mean stream 
velocity, small density ratio of solid and gas, etc.)“, the more 
uniform is the distribution of solid particles. 

The gravity effect has been neglected in the above considera- 
tion. In some of the results, however, its effect is quite obvious. 
Assymmetry of concentration profiles actually indicates super- 
position of the gravitational acceleration on the acceleration 
due to transport of momentum from the wall. The greater the 
value of f when compared to the gravitation of acceleration, 
the more symmetrical is the concentration distribution (compare 
Figures 7(c) to 7(b)). 


Conclusions 


The above results suggest that in a two-phase turbulent 
stream of high average stream velocity, low loading (within 
the carrying capacity of the stream) and with small particles 
(hence small gravity effect"), the concentration profile of 
solid particles in steady fully developed turbulence would ap- 
proximate the turbulent velocity profile due to the similarity of 
these transport processes. 


The assumption of a constant concentration distribution 
appears to be valid for the conditions of low loading, high 
average stream velocity and small solid particle size. 


The above analytical consideration of mass transport in two- 
phase turbulent pipe flow is a simple way of visualizing the 
problem of concentration distribution in that case. The para- 
meter fR?/Du is a simple measure of the distribution; a smaller 
parameter means a more uniform distribution. 
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mass of solid per unit volume) 





C, = Static concentration or number density of solid particles 

C, = Concentration at the center of pipe 

D = Particle diffusivity 

yj = Acceleration due to transport of particles across the tur- 
bulent velocity field of the stream 

J, = Zero order Bessel’s function 

K = A normalizing constant 

N = Counts per second 

n = Rate of motion of solid particles per unit area 

R- = Pipe radius 

r = radial coordinate 

u = Mean stream velocity 

up, = Mean particle velocity 

a = +¥fR?/Du, parameter of concentration distribution 
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IAL SECTION 


The Corrosion Behavior of Aluminum 





in Natural Waters’ 


HUGH P. GODARD’ 


This paper reviews the behavior of aluminum in 
contact with natural waters. The influence of water 
composition and service factors on corrosivity to 
aluminum are discussed. Laboratory pitting test 
methods are described and data are presented to 
support the existence of a cube root (in time) rate 
curve for the pitting of aluminum in water. Informa- 
tion on the pitting of aluminum in 17 Canadian waters 
is given and mention is also made of the pitting of 
aluminum in seawater and oil field brines. A statis- 
tical treatment of pitting data is presented. Finally, 
measures for extending the life of aluminum equip- 
ment subject to pitting by water are suggested. 


hile the corrosion behavior of aluminum in contact with 
various waters has been studied as far back as 1920, 
the new metal has not been much used for handling natural 
waters until recent years. Some waters tend to pit aluminum 
and this may have discouraged its application. A 10 year re- 
search program in these laboratories on the pitting corrosion 
behavior of aluminum alloys has led to a better understanding 
of the mechanism and of the environmental factors that influence 
pitting and particularly of the corrosion rates. This together 
with the application of statistics to pitting data now enables the 
prediction of performance from laboratory test results and from 
early service experience. Further, it is becoming apparent that 
for many applications pitting can be tolerated providing perfora- 
tion does not occur and equipment is therefore being designed 
to contain the pitting for the desired life period. 
\ glossary of water and pitting terms used in this paper is 


appended. 


General Corrosion Behavior of Aluminum in Surface 
Waters 


In general, the corrosion behavior of aluminum in surface 
waters is good, but it can vary from excellent to poor depending 
on the composition of the water and the conditions of service. 
Conditions of service include such factors as water velocity, 
temperature and oxygen content, which are often as important 
as the actual composition of the water. 


If aluminum is corroded by water the attack takes the form 
of pitting. There is no general corrosion or gradual thinning 
as occurs with steel. The number of pits formed and their rate 
of penetration depend on the composition of the water and on 
the conditions of service. Pitting that would occur in stagnant 
water is often prevented completely if the water is in constant 
motion, either by flowing or by turbulence due to agitation® 
The rate of penetration of pits decreases rapidly, and if the 


1Manuscript received November 12, 1959; accepted June 15, 1960. 
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metal thickness is adequate, perforation will not occur for a 
very long time even though an appreciable pit depth may be 
reached after a short period of service. For example, although 

a piece of aluminum sheet 0.020 in. thick, pre-etched to render 
it active, will perforate i in Kingston tap water in only a week, 
a sheet only 0.125 in. thick is not expected to perforate in less 
than 25 years. In a long term experiment still under way the 
maximum pit after five years was only 46 mils. 


Influence of Water Composition 


In general, soft waters are the least aggressive, and most 
can be safely handled in aluminum:*. This is in contrast to 
black iron, galvanized i iron and copper to which soft waters are 
more aggressive than hard waters. 


The factors which are believed to influence both the 
occurrence and rate of pitting of aluminum in stationary waters 
are: 4,5,6) 

pH 

total hardness 
chloride 
oxygen 
copper 


The effects of the various compositional factors are inter- 
related, and it 1s not yet possible on the basis of the limited 
data available to w eigh the importance of each, or to predict 
from a table of water composition alone, the influence of the 
water on aluminum. The copper content is believed to be 
important, but as copper in most surface waters is less than 0.2 
p-p-m. it is rarely determined and is thus usually unavailable. 
The influence of copper is dependent on pH, and is less pro- 
nounced in alkaline waters’. More copper can be tolerated in 
a soft water than in a hard one. Traces of mercury which are 
present in a few instances in some ground waters and in indus- 
trially polluted streams have a strong detrimental influence on 
aluminum. Mercury also is not determined in conventional 
water analysis. 


Determination of the Pitting Tendency of a Water 
To Aluminum 


(1) General 


The pitting tendency of a given water to aluminum under 
stagnant conditions can be determined in the laboratory in a 
fairly short time with a five gallon sample of water. A pitting 
rate curve is determined which can be related to curves for 
other waters for which there ts service experience. \ reasonable 
prediction can then be made as to the corrosivity of such waters 
and the probable life of a piece of aluminum equipment with a 
given wall thickness in it. Velocity tests can also be made to 
establish the influence of steady flow. 
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Figure 1—Influence of temperature on pitting. 


(2) Laboratory Pitting Tests 


(a) Stationary water. he determination of pitting rate 
curves for aluminum in stationary water in the laboratory has 
been described by Aziz and Godard‘ Strings of 10 coupons 
are immersed for several time periods. Experience has shown 
that values taken at one week, one month and three months 
give an approximate pitting rate curve and a good idea of the 
aggressivity of a given water to aluminum. 

By rough estimate, in the eight years since this technique 
was devel: ped i in these laboratories, several hundred pitting rate 
curves have been determined. The early work was done on 
Kingston tap water in a 300 gallon tank. The water was 
adjusted to 25°C. and a low rate of input and overflow was 
on to avoid bacterial growth. The water changed 

ontinuously about once a day so that there was no appreciable 
water movement. 

It has been found that pitting tests can be run in smaller 
containers and common household plastic pails have been used 
to test waters brought to the laboratory from other locations. 
It has been shown that the same values are obtained in Kingston 
tap water whether in a pail or in the large tank. It has also 
been shown that the pitting of aluminum in laboratory tap water 
was the same as its behavoir in Lake Ontario water taken 
directly from the lake at a point remote from the city? 


(b) Moving water. |n work on the pitting of aluminum 
in Lake Ontario water, samples were put directly into the lake, 
near the shore. These did not pit at all. As ‘the water was 
fairly rough at the test location it was reasoned that perhaps 
the constant motion of the water had prevented the initiation 
of pitting. 


\ laboratory test was devised which showed that the number 
of pits and their depth after a given time was greatly reduced 
with increased rate of flow, and that at 8 f.p.m., a rather low 
rate of motion, no pitting occurred at all while the coupons in 
the adjacent stationary water pitted normally’. 

More recently a new apparatus to test aluminum in moving 
water has been developed that can be applied to volumes of 
water as small as two gallons. It is thus applicable to the five 
gallon water samples commonly brought to the laboratory. The 
test is being applied to lengths of tubing as well as to elongated 
sheet samples. 

The apparatus consists of a simple re-circulating system in 
which the water is drawn up out of an inert polyethylene 
container by means of a stainless steel gear pump equipped 
with Teflon plastic gears and pumped chrough a vertical stand 
of the tubing under study and then returned to the container. 
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Figure 2—Influence of temperature 
on individual pit current. 


The pump is driven by a variable speed drive (0 to 400 r.p.m.) 
which in turn is driven by an electric motor. All connections 
are made with Tygon plastic tubings secured with hose clamps. 
Two marks are placed five feet apart on the lower and upper 
short sections of connecting Tygon tubing, and at water veloc- 
ities of up to | foot per second the flow rate is determined by 
measuring with a stop-watch the time it takes water to pass 
between these two marks. At high rates of flow the velocity 
is determined by measuring the time it takes to discharge a 
known volume of water through tubing of known internal 
diameter. Using this apparatus and 1/2 inch diameter tubing, 
water velocities of 7 feet per second have been obtained. 

This work has shown that extruded aluminum is more 
susceptible to pitting than sheet surfaces. 

Field experience has shown that at much higher velocities, 
in the order of 20 f.p.s., turbulence occurs especially at fittings 
which cause pitting of aluminum. 


(c) Influence of temperature. It has been found that 
as the temperature rises the pitting probability increases and 
the pitting rate decreases (See Figure 1). The pitting of 25 
aluminum in Kingston tap water has been studied up to 70°C. 
using an individual pit current technique. For this purpose a 
pit is machined from a disc of metal and re-inserted in the hole 
using a plastic insert to insulate the pit electrically from the 
sheet. Electrical leads fitted to the back of the disc containing 
the pit, and to the sheet are connected to a microammeter. This 
enables measurement of the actual current flow from the pit. 
Some pits develop a maximum current at about 40°C. at which 
point the rate is five times that at room temperature (See 
Figure 2). Above 40°C. the pit current drops sharply and, 
at 70°C., becomes less than at room temperature. Other pits 
have no such maxima and the pit current decreases linearly over 
the whole temperature range. Thus, above 40°C., the higher 
the temperature the longer is the life. This is the reason for 
the recent use of aluminum in domestic hot water tanks. 


Shape of Pitting Rate Curve in Fresh Waters 
An examination of the laboratory pitting data suggested 
that the maximum pit depth (d) was proportional to the cube 
root of time (t)°®. Thus the rate of penetration may be ex- 
pressed by the equation: 
d = Kt"s 


Ic is believed that the constant (A) will be a function of the 


168 The Canadian Journal of Chemical Engineering, October, 1960 


(ay 


IN MILS 


DEPTH 


PIT 


MAX/MUM 


Figure 3 


alloy an 
moveme 
root of t 
to fit da 
on Alecat 

The 
time of 
the spec 
however! 
duce an 
from th 
The val 
will be 
such as 


will be | 
error W 


period ¢ 
The err 
pit dept 
parison 


Int 
The 
has cor 
from a 


Que., 1 


shown 



















































PIT OEPTH 


MAXIMUM 





TIME IN MONTHS (1) 


Figure 3—Pitting rate curve for Arvida 57S aluminum pipe. 


alloy and of the water conditions (comp ysItion, temperature, 
movement, etc.). A plot of maximum pit depth against the cube 
root of time gives a straight line. The equation has been found 
to fit data obtained on many alloys in Kingston tap water and 
on Alcan 3S sheet in some 24 other Canadian waters. 

The cube root equation assumes that pitting initiates at the 
time of immersion which 1s the case in laboratory tests where 
the specimens are etched to render the surface active. In service 
however, an induction period (t,) may occur, which will intro- 
duce an error in the cube root plot. Because time is measured 
from the start of service, the value ¢ will be in error by t — t.. 
The value for t'’* which more accurately, should be (t — t,)'/ 
will be distorted and a comparison of values at two periods 
such as 


d, — d, = K[(te — ¢,)' — (4 — é,)') 


will be incorrect to an extent depending on the value of t,. The 
error will increase as ft, becomes larger and will decrease as the 
period of observation (t) increases (t — t, will then approach f). 
The error can be avoided by plotting the cube of the maximum 
pit depth (d*) against ¢ since this is the same relation. A com- 
parison of data at two time periods then becomes: 


d.3 —_ d,3 = K [(te - i) = (ty ee t,)) 
= K (te _ ti) 


In this case ¢, disappears. 

The most convincing data in support of the cube root relation 
has come from the examination of sections taken at intervals 
from a 2000 ft. line of 4 inch Alcan 57S water pipe at Arvida, 
Que., now 13 years old. 
shown in Table 1. 


The maximum pit depths found are 


A cube plot of these data is shown in Figure 3. 
Taking the d* = Kt equation, and calculating the expected 
maximum pit depth for 13 years from the 3.3 and 6 year results 


TABLE 1 
Max. Pir Deprus IN ALCAN 57S WaTER LINE AT ARVIDA, QUE. 


Period of Service Max. Depth 
Years (mils) 
0.4 25 
3.3 31 
6.0 40 
13.0 53 
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Figure 4—Pitting rate curves for 3S alloy sheet in 
Canadian waters—the extremes. 


one obtains a value of 56 mils, in excellent agreement with a 
maximum depth of 53 mils which was actually found on a 60 
foot length of pipe. This is cited to emphasize the value of 
knowing the pitting rate curve and as an example of what can 
be done with pitting data. 


Typical Pitting Curves in Canadian Fresh Waters 


Figure 4 shows the pitting rate curves for 3S alloy in the 
least aggressive and the most aggressive of 17 Canadian fresh 
waters tested to date in this laboratory together with one from 
an English source. The least aggressive water was taken from 
the municipal supply at Shawinigan (South) Que., while the 
most aggressive water was found to be that from Mossbank, 
Sask. As all specimens were etched before immersion pitting 
started at the time of immersion. 

The partial composition of these waters as determined on 
the samples is shown in Table 2. 

To obtain a pitting rate parameter that can be visualized, 
the time in weeks to obtain a 40 mil pit in 3S alloy has been 
taken. On this basis the 17 fresh waters tested vary from 2.6 
weeks for Mossbank, Sask., to 953 weeks for Shawinigan 
South, Que. The data for the waters tested are given in Table 3. 

No simple correlations appear in these data. The pitting 
aggressiveness of these waters is not related simply to pH, 
total hardness or copper content or to any obvious combinations 
of these. In general, however, the softer waters tend to be at 
the top of the table while the harder waters tend to be at the 
bottom. If pitting aggressiveness to aluminum is to be related 
to water compositions more detailed water analyses will be 
required, but comparisons will be difficult owing to the wide 
variation in compositional factors in waters. 


Shape of the Pitting Rate Curve in Seawater 


Tests on many aluminum alloys immersed in the sea at 
three locations, (Esquimalt, B.C., Halifax, N.S. and Harbor 
Island, N.C.) for periods of up to tive vears have shown that 
in sea water also the rate of penetration of pits decreases with 
time. However, the number of samples per time point has been 
smaller and control over the tests much less, with the result 
that there has been considerable scatter in the results and an 
accurate determination of the rate curve has not been possible. 

In one test made in Halifax harbour, with the co-operation 
of Naval Research Fstablishment, Dartmouth, a pitting rate 


TABLE 2 


PARTIAL COMPOSITIONS OF WATERS FROM 
SHAWINIGAN (SOUTH), QUE., AND MOssBANK, SAsK 


Least Aggressive Most Aggressive 


Shawinigan, Que Mossbank, Sask 


pH 7.1 7.8 

total solids 140 p.p.m 1038 p.p.m 
total hardness 73 i 555 

copper 0.04 ss | 0.005 
chloride - 20 o 
sulphate - 417 
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TABLE 


3 


ParTIAL COMPOSITIONS AND PITTING DATA FOR SEVENTEEN FRESH WATERS 
(In increasing order of pitting corrosivity) 


an Wks. to Hardness Copper 
Order Fest No. 40 Mils. Location pH p.p.m. nes. 
1 8 953 Shawinigan South, Que. 1a 73 0.04 
2 7 453 Shawinigan, Que. 7.4 18 0.04 
3 12 207 Crofton, B.C. 6.7 27 0.021 
4 6 205 Hamilton Bay, Ont. i353 205 0.003 
5 6 Kingston, Ont. 7.9 160 0.005 
6 5 175 Credit Valley, Ont. re 0 0.028 
7 4 147 Columbia River, B.C. 1:3 72 0.017 
S 16 83 Canvon Meadows, Alta. 7.9 169 0.005 
9 14 46 Regal Golf Course, Calgary, Alta. 8.1 331 0.007 
10 1 25 N. Sask. River, Drayton Valley, Alta. 8.1 267 0.11 
11 13 23 Peterborough, Ont. 0 86 0.012 
12 15 17 R. G. May Golf Course, Calgary, Alta. 7.9 218 0.002 
13 22 8g Billingham Beck, England 8.7 443 0.011 
14 24 6 Jasper, Alta. 8.2 196 0.007 
15 2 6 Lethbridge, Alta. 7.9 228 0.017 
16 9 4.4 South Saskatchewan River, Sask. 16 206 0.04 
17 10 2.6 Mossbank, Sask. 7.8 555 0.005 
curve was established for Alcan 65S-76 sheet immersed in sea- TABLE 4 


water for periods up to 14 months. 
cube root plot is shown in Figure §. 


For etched coupons the 


Shape of Pitting Rate Curve in Oilfield Brines 


This is a subject of considerable current interest in this 
laboratory. Although only five brines have been tested to date, 


several others are now under investigation, while samples of 


still others are being obtained. 


It should be pointed out at once that oil field brines differ 
very considerably from sea water. Sea water has a uniform 
composition and while the salinity or salt content varies from 
one ocean to another and also with the season the ratio of salts 
is constant and is accurately known. An average solids content 
for sea water is about 40,000 p-p.m. Oil field brines may contain 
six or more times this amount of solids and the individual 
components vary widely. For example, the compositions of 
the five brines tested to date are given in Table +. 


Oilfield brines may contain hydrogen sulphide, in which case 
they are designated in which case they are 
termed “sweet”. The presence of H.S in brines has a very 
strong influence on their corrosiveness to steel. Its influence 
on the corrosivity of brines to aluminum is not yet known 
although H2S by itself is not aggressive to aluminum. 


“sour” , or not, 


Oilfield brines are believed to contain variable and sometimes 
appreciable amounts of heavy metals although analyses are 
lacking. It is expected that heavy metals, if present, would 
have marked influence on the pitting behavior of aluminum. 

All five brines tested to date tended to pit aluminum and 
decreasing rate curves were obtained. Two of the five showed 
rates of the order found in fresh water and seawater while the 
other three showed much higher pitting rates, and are considered 


too aggressive to be handled in aluminum. 


(da 
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Figure 5—Pitting rate curve for 65S aluminum in seawater. 
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PARTIAL COMPOSITIONS OF SEVERAL OILFIELD BRINES 


Series | Test Total solids Cl SO, . 
No. No. pH p.p.m. p.p.m. | p.p.m, HS 
1 3 om 155,700 absent 
2 11 6.6 127,000 | 68,600 | 9,500 | present 
3 17 1.3 15,000 8,300 4 | absent 
4 18 rae 16,000 7,700 | 1,800 | present 
5 23 8.9 94,800 50,900 | 4,000 | present 


By coincidence, the three more aggressive brines all contained 
hydrogen sulphide while the two less aggressive brines did not. 
The three aggressive brines also differed from the other two, 
and from all the fresh waters tested, in that no pitting developed 
for the first week even on etched coupons. 


Treatment of Pitting Data 
(1) Theory of Extreme Values 

A theory of extreme values was developed by Gumbel"! to 
predict probable weather and flood excesses from limited data. 


The Gaussian probability curve for samples of data with random 
distribution is generally known. 
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Figure 6—Extreme value plot of Arvida pipe data. 
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‘The “extreme value’ is the last point on the tail of the curve 
to the right. Gumbel showed that under certain very general 
conditions in a series of data the maxima themselves followed a 
distribution curve, according to the formula: 

N=e*” 

Analysis of a large volume of data involving the counting 
and me asuring of tens of thousands of pits has shown that the 
pitting of aluminum in fresh water follows Gumbel’s equation, 
and thus extreme value statistics can be applied to aluminum 
way data? 

\ special t type of extreme value probability graph paper has 
been devised to permit direct plotting and extrapolation of the 
results. Treatment of data can be illustrated by taking the ex- 
ample of one of the pitting strings of 10 coupons tested in this 
laboratory. The deepest pit on each coupon is found and 
measured. These measurements (7) are arranged in order of 
depth, then are numbered from one to ten and the probability 
of occurrence ($) of each determined from the equation: 


n 


= 41 





‘The maximum pit depth on each coupon is now plotted on 
the special paper with depth as abscissa and probability as 
ordinate. If a straight line can be drawn through the points, 
conformance of the data with the theory can be assumed. 
Obviously the more data available the better the chance of 
obtaining a good straight line and thus the maximum depths 
from 20 samples will be better than 10 and so on. 

The use of an extreme value probability plot is illustrated by 
the pitting data from a 60 foot length of Alcan 57S water pipe 
after 13 years’ service at Arvida, Que. This example was cited 
earlier to support the cube root pitting equation. The pipe was 
cut into 20 lengths three feet long and the deepest pit in each 
determined. After placing in order of i increasing depth, and 
determination of probability, the data were plotted on the 
extreme value graph paper (See Figure 6). The first thing to 
be noticed is that a good straight line can be drawn through the 
points. This proves that the data can be treated by extreme 
value statistics. It will be recalled that the deepest pit by actual 
measurement was 53 mils. However, the line indicates that the 
predicted maximum pit depth for the 60 feet, based on the other 
19 sections and as defined in the appendix is 57 mils. 

The right hand ordinate of the paper carries a scale desig- 
nated ‘“‘return period’. The “‘return period” for any given pit 
depth 1 is the number of observations that must be made to find 
a pit of that depth. Conversely for any given number of samples 
the predicted maximum pit depth can be determined. 

Now the pipe was 2000 feet long and the unit section length 
was 3 feet. There were thus 667 such lengths in the entire line, 
and for a return period of 667 the predicted maximum pit depth 
is 90 mils. In other words, if the entire 2000 feet of line had 
been examined this should be the deepest pit. 

Thus an extreme value plot will give several valuable pieces 
of information. 


TABLE 
ALUMINUM PIPELINE 





(1) If a straight line is obtained the data fit the theory and can 
be treated by it. 

(2) Data which lie well away from the line are spurious and 
can be discarded. 

The predicted maximum pit depth for the area examined 

can be calculated. 

(4) The data can be extrapolated to predict the maximum pit 
depth for any selected larger area of surface. 

Another recent application of the extreme value plotting 
technique has been to determine the degree of uniformity of 
pitting rate for 50 consecutive lots of Alcan 2S sheet produced 
in one plant over a period of one year using the maximum pit 
depth found in Kingston tap water on a string of 10 coupons 
per lot after three months’ immersion. Although the deepest of 
the maxima was double the shallowest maximum, an extreme 
value plot gave a straight line and thus it was concluded all the 
material was part of one population and there was no real 
difference among the lots. 

Eldredge has applied the method with success to the treat- 
ment of pitting data obtained by caliper surveys on steel oil 
well piping. 


(2) Calculation of Probable Life 

If the predicted maximum pit depth is calculated and the 
cube root pitting rate law assumed, the time to initial perforation 
can easily be computed. In the case of the Arvida pipe, 36 
years would be required to reach the pipe wall thickness of 125 
mils, or 23 more years’ service could be expected. From some- 
what limited experience the author’s present view is that there 
is a factor of safety in such predictions, and that an even longer 
life will probably be realized. 


Case Histories 


(1) Water Pipelines 
The data in Table 5 were compiled during regular examina- 
tions of actual aluminum pipelines. 


(2) Water Storage Tanks 
At Shawinigan, Que., two large (50,000 gallon) aluminum 
overhead water storage tanks were built in 1948, one for an 
industrial plant, and one for the municipality. When inspected 
in 1957 after 9 years’ service the maximum pit in the first 
mentioned was found to be 90 mils and in the other 184 mils. 
Calculations of probable life predicts a long forward life free of 
maintenance at least 100 years for the industrial plant tank, 
and 25 years for the other. 


(3) Heat Exchangers 

(a) Wall board factory installation. |n a large mill 
near Ottawa an all-aluminum heat exchanger is used to recover 
waste heat expelled as steam during the pressing of a waste- 
wood type of wall board. The steam is drawn by a fan into a 
hood and then through a bank of water sprays (Ottawa River 
water) to warm up the i incoming plant water supply. The spray 
chamber is about + X 5 X 10 feet, and there is a series of baffles 
to remove entrained water droplets. The exchanger was placed 


5 


Case HisTories 


(ae pit depths found in service) 





Picccsteil —— Years 
escription -oca tion Service 
a sensi nntisteatinta = 
Plant Water Supply Arvida, Que. 13 
Golf Course Supply | Kingston, Ont. 11 
Golf Course Supply | Calgary, Alta. 7 
Golf Course Supply | Peterborough, Ont. 6 
Goif Course Supply | Joliette, Que. 5 
Steam Condensate Return | Kingston, Ont. 12 
Water Injection Line Drayton Valley, Alta. | 1 


oe i a Pipe Maximum Length 
Length | Diameter | Wall | Pit Depth | Examined 
= — Mils Mils Feet 
2000 | 125 53 60 
300 114 to2% 125 0 60 
2 55 37 40 
3 55 15 20 
2,3,4 57 Perforation 10 
500 1% 133 0 24 
39,500 4,6 188 andup 6 3 
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. bafllles removes entrained spray. 


in service in 1949 to replace a similar unit of galvanized steel 
that had lasted only two years. When examined one month 
later it was in excellent condition, and it is-still in use after 9 
years service without maintenance. 


(b) Air conditioning system. In a large Canadian 
manufacturing plant an all-aluminum heat exchanger is used to 
de-humidify incoming air to the main air-conditioning system. 
The air is drawn through a bank of water sprays by a large fan. 
The spray chamber 1s ‘about 12 X 12 X 12 feet. A series of 
The unit was placed in service 
in 1950. When examined in 1951 it was in excellent condition. 
It is still in use after 8 years service. 

The only corrosion problem encountered was perforation of 
a horizontal sheet in the fan plenum where water which had not 
been removed by the baffles settled in stagnant pools. The 
violent agitation in the main chamber apparently prevented such 
pitting. 


Prevention of Perforation of Aluminum Equipment 


The following means are available for extending the life of 
aluminum equipment handling a water which tends to cause 
pitting: (1) Increase the wall thickness. The value of increased 
wall thickness is apparent from the nature of the pitting rate 
curve. Twice the wall thickness will give eight times the life. 
(2) Use Alclad aluminum. This consists of a duplex material 
the core alloy of which is protected by a thin cladding layer of 
another aluminum alloy selected to be anodic to the core so that 
when a pit reaches the core it will spread laterally. Further 
penetration will be arrested until the area of core exposed is 
such that the cladding will no longer protect the center due to 
the increased resistance of the electrolyte patch. Experience 
indicates that this delaying action will increase the time to 
perforation of a given thickness of metal from five to ten times. 
(3) Increase the movement of the water, i.e. avoid stagnant 
conditions. This has already been dealt with. (4) Deaerate the 
water. It has been shown that deaeration of Kingston tap water 
prevents the initiation of pitting that would otherwise occur 
and it seems reasonable to assume that deaeration will have a 
similar influence in other waters. (5) Use inhibitors. In a 
circulating system an inhibitor such as sodium dichromate at 
500 p.p.m. and pH 8 should completely prevent pitting. There 
is some evidence that poly phosphate and silicate inhibitors also 
inhibit the pitting of aluminum in water. (6) Apply cathodic 
protection. It has been shown that the pitting of aluminum in 
Kingston tap water can be prevented by the initial application 
of a protective current density of about 2 Ma/ft* at a potential of 
about 1.0 volts negative to a copper-copper sulphate electrode, 
which is the value necessary for the prevention of pitting of 
aluminum in seawater. (7) Avoid undesirable dissimilar metal 
couples. Aluminum in contact with dissimilar metals other 
than zinc and cadmium tends to corrode galvanically, the severity 
depending on the conductivity of the water, and on the other 
metal. In fresh waters the tendency is small but in brines and 
seawater it is large. The most aggressive metals are copper, 
brass and bronze, followed by lead and then steel. Stainless 
steel can be safely used with aluminum in all but brines and 
seawater. 


Conclusion 


It has been shown that while some waters tend to pit 
aluminum, the rate of penetration follows a rapidly decreasing 
rate curve that approximates a cube root function. In many 
cases such pitting can be tolerated since in water handling 
equipment perforation i is the main criterion of performance. It 
has also been shown that statistical treatment of pitting data is 
possible, and is, in fact, a valuable tool in predicting the life of 
equipment from early service experience. The amount of 
aluminum equipment in water handling service is increasing 
rapidly and should continue to grow as the performance of 
aluminum becomes more widely known. 
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Nomenclature 


d = maximum pit depth 

constant 

time of service 

time to initiation of pitting 

maximum pit depth at time ¢, 

maximum pit depth at time fy 

probability 

ordered number of the pit depth (1 to 10) 
Total number of maxima (10) 


PTS QAacr yy 
tit an ea 


= 
I 
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APPENDIX 


Glossary of Terms 


(1) Water Technology 


Natural waters. Vhe term “natural water’ refers to 
rain water, surface water and ground water. These are usually 
fresh waters, as opposed to seawater. The term “surface water” 
describes waters in creeks, rivers, ponds and lakes, while the 
term “ground waters” denotes sub-surface waters from springs 
and wells, some of which are brines. Any of these waters 
which have been processed for consumption or industrial use 
are termed “treated waters” 

The ~ water found in nature is dew, hoarfrost and 
rainwater. However, many rain waters are far from pure 
because rain drops falling through the atmosphere absorb matter 
in amounts depending on the purity of the air. Above large 
manufacturing areas the air is often laden with dust, soot, silica, 
silicates, sulphates, chlorides and carbonates, oxides of sulphur 
and nitrogen, hydrogen sulphide, ammonia, organic matter and 
other impurities. T hese are absorbed in the falling rain or snow 
and give rise to very inferior water. Near the seacoast, rain 
carries salt absorbed from particles of salt or sea spray blown 
inland by on-shore breezes. This may reach 15 p.p.m. but 
decreases rapidly with distance from the sea. 

Surface waters dissolve salts from the soils with which they 
are in contact and thus have varying composition depending on 
the soil and the length of contact. 


Water composition units. No standard method of 
reporting the results of water analyses has been adoped in 
Canada or elsewhere. Several systems are in vogue necessitat- 
ing the use of conversion factors when comparing analyses. 

Constituents may be reported as grains per gallon, (Imperial 
or American) in parts per hundred thousand, or in parts per 
million. In Europe results are sometimes reported as milligrams 
per litre (identical with parts per million). 

Some laboratories report the bases as basic radicals, others 
as oxides, and the acids may be reported similarly or as hypos 
thetical combinations. The system used by the United State- 
Geological Service is widely used, has been adopted by the 
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parts pc! million. The following items are commonly reported: 


pH iron 
residue On evaporation silica 

total hardness bicarbonate 
calcium sulphate 
magnesium nitrate 


alkalis (sodium and potassium) chloride 

Suspended matter, turbidity and color are sometimes re- 
corded. Analyses may also be made for dissolved oxygen and 
free carbon dioxide, but these require determination at the time 


of sampling. 


Hardness in fresh water. Originally “hardness” was 
defined as the capacity of a water for precipitating soap. 
The fatty esters such as sodium oleate, stearate, etc., instead 
of remaining in solution and forming a lather are precipitated by 
calcium and magnesium (present as the bicarbonates, sulphates 
and chlorides) as an insoluble curd. Other salts such as those 
of iron, aluminum, manganese, strontium, barium, heavy metals, 
brines and free acids precipitate soap in a similar manner. 
Normally, these salts are present in insignificant amounts and 
the hardness reported is due only to the calcium and magnesium 
salts. This is known as “carbonate hardness” or “temporary 
hardness”. The hardness due to sulphates and chlorides is 
known as “‘non-carbonate”’ or “permanent” hardness. 

‘Today, hardness can be titrated quickly with a visual color 
change end-point using reagents based on EDTA (ethylene 
dianane tetra acetate). In this method small amounts of the 
salts producing permanent hardness are reported as calcium, an 
insignificant error in most cases. 

The simplest general compositional classification of fresh 
waters is by hardness, which is expressed in terms of calcium 
carbonate in parts per million. Waters are classed as soft or 
hard according to the following table: 


very soft 1— 30 p.p.m. (CaCO) 
soft 31— 60 
medium hard 61-120 
hard 121-180 


very hard > 181 
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Canadian Department of Mines and Technical Surveys, and is 
used at Aluminium Laboratories. The unit in this system is 


Other hardness units, and their conversion factors to p.p.m. 
are given in the table below: 


X Factor 
Unit to p.p.m. 
parts per 100,000 10.0 
grains per U.S. gallon 17.1 
grains per Imp. gallon 14.3 
English (Clark) Degrees 14.3 
French degrees - 10.0 
German degrees 17.9 
Equivalents per million 50.0 
nglish Degree = 1 grain CaCO; per Imp. Gal. 


= 10 p.p.m. CaCO; 


1 En 
1 French Degree 
1G = 10 p.p.m. CaO 


rerman Degree 


(2) Pitting Terminology 
Pitting rate. The rate of deepening of a pit, which 
follows a decreasing rate curve in the case of aluminum. 


Maximum pit depth. The depth of the deepest pit 
found on the surface examined. Not necessarily the deepest 
: : : P 

pit on the entire surface of the equipment. 


Predicted maximum pit depth. he most probable 
value for the deepest pit on the area being considered, predicted 
by extreme value statistics from actual pitting data from the 
same or similar samples. 


Pitting probability. he tendency to initiate pitting. 
May be expressed as a percentage from test data‘:?). 


Pitting density. Vhe number of pits per unit area. 


Severity of pitting; Intensity of pitting. Ambiguous 
terms which- should not be employed since it is not clear 
whether they refer to pitting density or pit depth. 


Average pit depth. Frequently the five or ten deepest 
pit depths are measured and averaged. Such numbers have 
little practical value since only the deepest pit is of concern. 
Much less useful than predicted maximum pit depth. 


x * * 






























